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Final Report on Delivery Order 48 of NAS8-36955 

This contract supported four activities. 

1 . Differential Collision Cross-Sections 

The primary purpose of this effort was to provide on-going support for Dr. 
Charles Keffer to establish a facility at MSFC for the measurement of differential 
collision cross-sections of O on N 2 and other gases of relevance to understanding 
vehicle-environmental contamination effects in orbit. 

This work was initiated in 1988 under NAG8-100, entitled "Space Vehicle 
Contamination Study." During the period 6/1/89 to 11/6/90 the work was funded 
under D.O. 48, D.O. 59 and NAS8-37586. In June 1990, Dr. Keffer was awarded a 
grant with himself as Principal Investigator and has continued the work since that time 
in his own capacity without further assistance from D. G. Torr. 

To place the flow of funding in perspective, Table 1 shows the percentage of Dr. 
Keffer’s time charged to the above mentioned contracts in support of the fabrication of 
the cross-section facility. Also to provide some traceability of progress, we review 
briefly below the overall chronological evolution of the work, identifying what was 
specifically accomplished under D.O. 48. Also Dr. Keffer provided weekly written 
reports on his progress consistently and these are available on file. 

The proposed experiments use a 5 eV atomic oxygen beam for differential 
scattering cross-section measurements. This 5 eV O beam is formed by focussing the 
output from a high energy pulsed C0 2 laser into the throat of a molecular oxygen nozzle 
beam source. 

NAG8-1 00 

Work on the cross-section facility was initiated under this grant in 1988. The 
final report on progress made was submitted on November 7, 1989. Essentially the 
conceptual design for the facility was developed, machining and welding were completed, 
and documentation generated. Signal to noise ratio estimates were calculated and writing 
of software for operation of the instrument initiated. 

N AS8-37586 

This contract supported Dr. Keffer from June 1 through September 30, 1989. 
During this period the cross-section chamber was assembled, and initial problems with 
the system debugged. By the end of the 4 month funding period the system was pumped 
down, but it was discovered that the chamber went out of alignment under vacuum 
conditions. Work on rectifying this problem was initiated under D.O. 48. 

Work Accomplished Under D.O. 48 

This phase started with an evaluation of the noise sources in the system, and in 
correcting the serious alignment problems identified under NAS8-37586. 

The correction approach that was adopted was to adjust the assemblies attached to 
the bellow section. Following the alignment correction, the chamber was successfully 
aligned under vacuum on one axis. Vent and evacuation cycles were performed for 
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alignment reproducibility. The alignment of the CO 2 beam was successful to the point of 
hitting the beam valve aperture with the focussed laser beam. However, the coating for 
the mirror which reflects the CO 2 laser beam was destroyed because of the high energy 
density of the beam. 

Various materials were tested for CO 2 laser mirrors, and nickel copper mirrors 
with a gold coating were found to be successful. 

Alignment tests were performed on the mass spectrometer housing, and work 
commenced with thermal scattering measurements. The testing indicated that all critical 
elements could be aligned. Estimates were made of the accuracy and uncertainty of the 
alignment. This was the status of the project on January 31, 1990. 

During February work continued under D.O. 59. Under D.O. 59 work began on 
detecting a scattered signal from O 2 on N 2 collisions. However, the signal to noise was 
poor. Shielding of signal lines corrected the problem. 

Work on the system operating software was re-initiated. 

Work was started in a pulse mode of operation with crossed beams. 

An anomalous result was obtained for the O 2 -N 2 beams. 

At this point funding was continued under D.O. 48. 

The problem causing the anomalous results was identified and corrected and both 
beams exhibited sharp temporal and spatial profiles. 

To increase sensitivity further the filament emission of the mass spectrometer 
was increased to give a factor of 10 improvement in detection efficiency. Also, the 
detector aperture was increased. 

An anomaly in the results obtained was discovered in early March. This was 
eventually traced to a grounding problem in the instrument triggering circuit. 

Action taken to improve S/N included: 

Inhibition of the two pulsed valves controlling the sources of the crossed beams 
when data is not being taken. 

An RF shield for the mass spectrometer data cable was added. 

This was the status of the project at the end of April 1990, after which the 
project was funded under D.O. 59 until Dr. Keffer obtained independent funding to 
continue the work. 

2 . Groundbased Scientific Observations of Rocket Releases during 

NICARE-1 

This project utilized the broadband imaging spectrometer referred to as the 
Compact Spectrometer developed in a colloborative effort with MSFC under NAG8-060. 
Since the Compact Spectrometer was the only instrument in existence capable of 
measuring a spectrum simultaneously across the wavelength range 3600 - 11,000 A, 
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UAH was approached by the Naval Research Labs (NRL) to monitor a rocket launched 
chemical release experiment in which the emissions of two new chemicals were to be 
identified. The spectra to be expected were essentially unknown. Thus the need for a 
broadband spectral imager. 

The measurements were made from a naval station at Duck, North Carolina on 
October 22, 1989. The rocket was launched from Wallops Island, Virginia. Figures la 
show the twilight airglow spectrum recorded prior to launch. Figure 1 is the nightglow 
spectrum recorded two minutes prior to launch, and Figure 2, the same one minute after 
the flight. Figure 3 shows the spectral image recorded for 30 seconds immediately prior 
to the first release, which was a nickel carbonyl release. The natural airglow 
background has been removed. Figure 4 shows the spectral image recorded for 30 
seconds commencing after the nickel carbonyl release. Figure 5 shows the same 60 
seconds after the release. Emissions were only observed in the field of view of the mid- 
visible channel indicating that a sharp edge must characterize the cloud. The emissions 
appear to be a CO 2 continuum. Figures 6 and 7 respectively show the corresponding 
spatial distribution of the CO 2 continuum. The narrow width of the emission suggests 
that the rocket exhaust plume was being observed, and that no emissions from the 
release were observed. No emission was detected from release 2 which was a 
trifluromethyl bromide release. Figures 8 and 9 show a blow-up of the two CO 2 spectra. 

These data were provided to NRL and the results were presented by Dr. Paul 
Bernhardt at the Fall 1 990 Meeting of the American Geophysical Union. A copy of the 
printed abstract is attached. 

3 . Data Compression Study for the UVI 

The purpose of this study was to evaluate ways by which the data rate generated 
by the detector readout could be matched to the telemetry transfer rate to the ground. 
The main goal was to determine what constraints the science requirements placed on 
temporal and spatial resolution, because the data rate can be reduced by increasing the 
CCD on chip integration time and by reducing spatial resolution requirements. Also 
placement (mapping) of the object dimensions in the image plane, i.e. a round earth on a 
rectangular image plane provides some additional freedom. 

The pixel dimensions of the CCD are 244 x 550 = 1.342 x 10 5 pixels per 
readout, digitized to 12 bits yields a data rate of 1.5 x 10 6 bits per readout. If no spatial 
compression could be tolerated, this would dictate an on-chip integration time of -134 
seconds to match the 12 kbps telemetry rate assigned to the UVI. 

An evaluation of auroral dynamical time scales revealed two boundaries in this 
temporal regime: less than 1 second and greater than several tens of seconds. A thirty 
second integration time therefore requires a factor of 4.5 additional compression. The 
desired spatial resolution should be better than 40 km. 

A factor of 2 can be obtained by noting that the pixel size of the Thomson TH 7866 
CCD is 27 x 16 microns, which yields an undesirable elongated rectangular spatial 
resolution element. Thus it was decided to sum two of the 16 pm pixels to give a spatial 
resolution element of 32 x 27 pm as indicated in Figure 3.1. This image is fiber 
optically mapped from the CCD to the focal plane of the image intensifier as illustrated in 
Figure 3.2, yielding a focal plane spatial resolution element of 87.4 x 73.7 pm which 
matches the optical blur spot size of the instrument. The ground spatial resolution 
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achieved at 9 Earth radii is 25 km which meets the science requirements. Thus an 
additional factor of 2.2 compression is needed to match the telemetry rate. 

This can be readily achieved with standard data compression techniques. For 
example, the first pixel in each row can be digitized to 12 bits, and then the intensities 
recorded with respect to the first expressed on a relative scale at lower quantization. 
Eleven bit digitization, for example yields a compression of a factor of 2. As a check for 
self-consistency the last pixel will also be digitized to 12 bits. 

This study resolved the data compression issue. Figure 3.3 illustrates the 
placement of the image on the image tube focal plane. 

4. Science Priorities for Ultraviolet Imaging in the Mid-1990's 

This task was approached mainly from the perspective of the ionospheric 
priorities for global imaging, and not from the magnetospheric perspective because our 
expertise lies in ionospheric physics. The magnetospheric input was provided by other 
Co-investigators with expertise in magnetospheric physics. 

The Earth's upper atmosphere is constantly undergoing periodic and aperiodic 
changes induced by solar eletromagnetic and corpuscular radiation. 

The primary energy source is the absorption of solar ultraviolet radiation which 
under quiet conditions produces predictable cyclic patterns of behavior. 

The solar wind introduces plasma into the terrestrial environment which 
energizes and penetrates into the Earth's magnetosphere where it is stored and then 
impulsively released creating, for example, the visual aurora. 

These processes generate complex and variable electric potentials and currents 
which, in turn, result in large insertions of electrical and corpuscular energy to the 
upper atmosphere, greatly perturbing is quiet time behavior. 

Figure 4.1 illustrates the energy sources to the ionosphere/thermosphere. Our 
goal is to determine the energy sources and their impact on the system. 

To understand the significance of the role of the UVI Imager for the ISTP mission, 
one can pose the questions: what will the UVI do that has never been done before. The UVI 
will provide global coherent images of terrestrial emissions every 30 seconds, i.e. a full 
global snapshot of the entire earth repeated at 7 wavelengths. This has never been 
achieved before, because previous images of the earth had to be built up one pixel at at 
time by scanning a pencil beam field-of-view across the earth. It required 15 minutes 
on the Dynamics Explorer to build up single image at one color. Thus the coherence 
needed over time scales of the order of tens of seconds to understand the global behavior 
of the dynamic aurora could not be realized. 

Second, the filter technology did not exist prior to UVI to separate out the 
ultraviolet wavelengths or colors emitted by the upper atmosphere, so that the full 
information of the spectral nature of the emissions could not be realized. (Note the 
dayside of the earth can only be viewed at VUV wavelengths). Thus the questions reduced 
to: what will coherent spectrally pure global imaging achieve? 

To answer this question we examined the status of ionospheric research. 
Currently models exist which have the potential to describe the behavior of the global 
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ionosphere. These ionospheric global models require the following fundamental input 
parameters: 

Solar EUV flux 

Global precipitated particle flux 

Upper atmospheric neutral winds 

Global electric field pattern 

Knowledge of the neutral atmosphere, i.e. composition and temperature 

If the ionospheric models are combined with a thermospheric global circulation 
model (TGCM) then the required inputs reduce to 

Solar EUV flux 

Global precipitated particle flux 

Global electric field pattern 

Under this contract we developed procedures to demonstrate that the UVI images can be 
used to provide the global precipitated particle flux and the global electric field pattern. 
These two parameters have been impossible to acquire on a global grid in the past, 
because the data are only useful if near simultaneous measurements can be made on a 
global scale. 

Under D.O. 48 we demonstrated that coherent global imaging of spectrally pure 
emissions will provide the data needed to determine the particle fluxes and electric 
fields. 


The first task to be performed was to determine what emissions could be used to 
characterize the particle energy influx. Using the auroral deposition code developed 
under NAS8-37586, the response of several candidate emissions to auroral particle 
energy input was evaluated. It was demonstrated that the molecular nitrogen Lyman- 
Birge-Hopfield bands in the Vacuum Ultraviolet (VUV) range 1200 - 1800 A (termed 
LBH U ong) vary linearly with the total precipitated particle energy flux. No other 
emissions exhibit this behavior. The results of these calculations are shown in Figure 
4.2. These results were presented to the UVI Science Team on August 15-16, 1990. 
This success was fundamental to the science objectives that could be developed, and can 
be attributed directly to the successful development of a rectangular VUV filter under 
NAS8-37586 that straddled the LBH LO ng wavelength range. A full scientific under- 
standing was developed as to why the LBH U ong wavelengths respond in the desired way. 

The next discovery proved to be equally significant, namely the ratio of two 
emissions were found which were sensitive only to the energy spectrum of the incoming 
particles, as defined by a "characteristic energy". Figure 4.3 shows the characteristic 
energy as a function of the ratio of the LBH LO ng/LBH S hort emissions, where LBH short lie 
in the wavelength range 1400-1600 A . Once again, the key to achieving this success is 
the ability to fabricate a filter that can spectrally isolate the LBH LO ng and LBH SHO rt 
emissions. The ratios of many emissions were investigated, and in each case with the 
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exception of LBH LO ng/LBH S hort. strong sensitivity to atmospheric composition and 
temperature were observed, which impacted their usefulness. 

The above findings mean that the spectrum and energy flux of the incoming 
particles will be derived from the UVI images on a global scale, every 30 seconds. The 
impact of this data on modeling will be profound for the following reasons. The auroral 
code, referred to above, can now be used to calculate the global ionospheric corpuscular 
source function as a function of time. This will be the first time that real values for this 
critical ionospheric input parameter will be available. The corpuscular source accounts 
for about half the total energy input to the ionosphere and thermosphere and in the 
auroral zone it is the dominant energy source. 

The results of this work were published in the Journal of Geophysical Research 
under the title, "The Dependence of 01 1356 and N 2 LBH Auroral Emissions on the 
Neutral Atmosphere", which is attached as Appendix A. 

Figure 4.1 illustrates the energy input sources to the thermosphere/ionosphere 
system. An important goal in Solar Terrestrial Physics is to quantify the effects of 
coupling of magnetospheric electrical energy into the ionosphere. Under this contract 
we developed a strategy to use the imager to help define the global electric field system. 

Current systems in the magnetosphere generated by interactions with the solar 
wind generate a cross-tail potential which maps down into the ionosphere, where it 
accelerates the ions. The ions collisionally impart momentum to the neutral atmosphere 
which sets up winds which typically can reach velocities of the order of kilometers per 
second. Thermally driven winds can also act on the ions generating ion drift motions. 
The effects of electric field are therefore manifested in ion and neutral motion. The 
former can be measured by groundbased incoherent scatter radars and the latter by 
Fabry-Perot interferometers which monitor a bright airglow feature at 6300 A. Also, 
the current systems associated with the electric fields cause fluctuations in the earth’s 
surface magnetic fields which can be observed by magnetometers. In addition, the 
electric field instruments on the ISTP spacecraft make direct in situ measurements of 
electric fields. All the above techniques can provide valuable spot (isolated) 
measurements of the global electric field. However, the imaging data can also be used to 
infer electric fields on a much broader scale. This statement is based on the fact that the 
ionosphere constitutes an important component of the magnetospheric electrical 
circuits. Currents can flow into and out of the ionosphere if the input and output points 
are electrically connected horizontally. The current carriers that close the loop are the 
electron-ion pairs that constitute the ionospheric plasma. Thus when auroral electron 
precipitation occurs, generating the auroral emissions, there is an associated production 
of electron-ion pairs. This production rate can be computed by the auroral deposition 
code from the UVI images as described above. Naturally, at night, data taken by the Iowa 
Visible Imager can also be used. This source function coupled with an ionospheric global 
model can be used to calculate the ionospheric electron density, i.e. the density of the 
current carriers, which means that the electrical conductivity of the circuit can be 
computed. These conductivities can then be used in an ionospheric electrodynamics 
model to infer the input electric fields. 

Currently large scale potential patterns describing ionospheric convection 
patterns are estimated using statistical models for the ionospheric electrodynamics. One 
such model referred to as the Kamide-Richmond model, utilizes the following input data 
to infer the large scale potential pattern: 
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Incoherent scatter radar convection observations 
Satellite direct field observations 
Ground and satellite magnetometer observations 
Ionosphere electrical conductance observations 

The groundbased observations are capable of tracking rapid changes in electric 
fields and currents, conductivities and associated magnetic perturbations. The 
electrodynamics models can combine different kinds of electrodynamical data to infer the 
high-latitude electric potential patterns as they evolve in time. However , spatial 
coverage is limited so that many stations would be needed to obtain global coverage. By 
using the global conductances derived from the UVI/VIS images the input database can be 
greatly enhanced, and at much higher time resolution, allowing the time-evolution of the 
global ionosphere to be modeled. 

Figures 4.4 - 4.9 illustrate our strategy for reaching a quantitative global 
modeling capability by the mid-1990's. 

Figure 4.4. shows the basic starting point, discussed above. P. G. Richards and D. 
G. Torr have developed a global ionospheric model referred to as the Field Line 
Interhemispheric Plasma (FLIP) model. This model computes the global electron 
density as a function of height and time. It requires several inputs which are identified 
in Figure 4.5. As is evident from this figure, the FLIP model produces a large number of 
outputs. Figure 4.6 includes the need for these key inputs in the flow chart. 

In addition to the auroral source function, the FLIP model requires as input the 
solar EUV flux, neutral winds and the neutral atmosphere. There are several ways to 
obtain these. The first option is to use empirical models which provide all three inputs 
as a function of geophysical indices such as magnetic and solar activity. However, these 
are essentially statistical models derived from existing databases, and do not necessarily 
accurately represent conditions on any specific day. The empirical models can be 
further constrained by groundbased airglow measurements as illustrated in Figure 4.7. 
UAH (our team) has been deeply involved over the last five years in developing 
techniques for deriving the solar EUV flux, the neutral atmosphere and winds from 
groundbased measurements of airglow and electron densities. 

To include the effects of electric fields, the FLIP model is coupled to an 
ionospheric electrodynamics model (see Figure 4.8), which derives input from the 
groundbased network of radars, ionosondes, magnetometer chains and in situ ISTP E 
fields. The electrodynamics model provides the electric fields to the FLIP model which 
returns the global conductances to the electrodynamics model, thereby iteratively 
improving the global electric field calculation. The FLIP model predicts the electron 
density profile as a function of time at the locations of ~ 100 ionosonde and ~ 6 to 7 
radar stations. Iteration is continued until the best electron density fit is obtained at all 
stations. 

Finally, the entire system is coupled with a global thermospheric circulation 
model (TGCM) as illustrated in Figure 4.9. The TGCM utilizes the joule heat source and 
ion drag momentum source from the ionospheric models, and provides the neutral 
thermosphere and winds back to the ionospheric models. Iteration is continued until the 
global electron densities best fit the groundbased measurements. 
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We believe that the science priorities developed here will significantly advance 
our understanding of thermospheric, ionospheric magnetospheric coupling. The role of 
the ISTP database will be to provide near-continuous global coverage of key inputs to the 
system on a 30 second temporal grid allowing the actual response of the high and low 
latitude ionosphere and thermosphere to all significant energy sources to modeled in a 
time evolution manner. 

This scheme is a very ambitious one, requiring large computing resources and 
excellent coordination between NASA and the NSF's groundbased network of aeronomy 
stations. However, it is highly viable. It can be done, and should be done. 

Reprints of relevant papers published on this work are attached in Appendix A. 

5 . Assessment of Optimizations Possible in UV Imaging Systems 

In this study, NASA's future needs in ultraviolet imaging were assessed, and 
possible optimizations evaluated. The following instrumental requirements were 
identified for the post ISTP area in order to meet the scientific goals of the future. 

1 . A capability will be needed to image the weak aurora from space. 

2. It was noted that strategic planning activities by NASA Headquarters 
during the late 90's had identified imaging of magnetospheric extreme ultraviolet 
emissions as a key goal for the future. At that time, however, the technology did not 
exist to provide the spectral purity and sensitivity needed to successfully image the key 
ions of the magnetosphere. Specifically, emissions of ionized helium and atomic oxygen 
at 304 A and 834 A respectively were identified as the prime emissions for imaging. 
Two-dimensional imaging of the magnetosphere could provide the first coherent 3-D 
picture of the ion populations by tomographic inversion of the images. 

The above two goals translate into the following new technology requirements: 
Better VUV filters than those developed for ISTP 

A new class of instrument is needed with 100 times greater sensitivity 
than the existing EUV capability, with 50 A spectral filtering. 

A study was carried out to investigate ways of increasing further the 
performance achieved with the ISTP filters developed mainly under NAS8-37586. This 
problem was solved with two new technology developments: 

1 . The development of n multilayer coatings - a new design concept. 

2. The use of multiple reflectors to improve resolution. 

A detailed report on the approach has been submitted for publication in Applied 
Optics under the title, "Multiple Reflectors as Narrowband and Broadband VUV Filters". 

This technology resulted in a decrease in bandwidth of the ISTP filter at 1356 A, 
for example, from ~5.0 nm to 3.2 nm with an increase in transmission from ~25% to 
~40%. The work was also presented at the Fall Meeting of the American Geophysical 
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Union in San Francisco, California, December 1990. A copy of these reports is attached 
as Appendix B. 

In addressing the question of extreme ultraviolet imaging, the thin film design 
program developed for the VUV was applied to the EUV region at 83.4 nm. Materials 
were identified which could serve as multilayer coatings in this wavelength regime. 

Coatings with 60% reflectance at 83.4 nm with out-of-band reflectance of 5% 
were designed. These coatings were incorporated into an EUV imager design in a highly 
innovative way; the design requires no actual filters in the system. The ISTP UVI 
instrument was selected for the imager optical design. This is a three-mirror system. 
It was found that by placing the reflective filter coating directly on the mirror surfaces, 
the mirrors could be used simultaneously as both the imaging and filtering elements. 
The design yielded a net system efficiency of -20% compared with -0.2% achieved with 
existing technology, thus providing two orders of magnitude improvement in sensitivity, 
which meets the sensitivity specifications for magnetospheric imaging of 0 + ions. 

This work was reported in a CSPAR presentation to the University President, Mr. 
J. Moquin, on January 22, 1991. Copies of the relevant charts used are included in 
Appendix C. 
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' SA51A-8 1052h 

Evolution of Three-Dimensional Plasma Cloud* in the 
Ionosphere 

T-Z Ma and R W Schunk (Center for Atmospheric and Space 
Science, Utah State University, Logan, Utah 84322-4405) 

A three-dimensional time-dependent model for the motion of 
plasma clouds in the ionosphere has been developed. In the 
model, the ions created from the released neutral gas and the 
background ions are treated separately so that the accelera- 
tion of the background ions and the deceleration of the cloud 
ions can be calculated. The model takes account or the finite 
conductivity along the magnetic field, pressure gradients, stress 
tensor, inertiakeffects, collisions, thermospheric winds, and cur- 
rent coupling between the cloud and background ionosphere. 
Studies are being conducted to contrast the expansion char- 
acteristics of both negative and positive ion clouds, long thin 
cloudlets and nearly spherical clouds, clouds with and without 
background neutral winds, and clouds that are both preionized 
and photoionized . Some of the cloud expansion features will be 
presented. 


SA51A-9 1106h 

N umeral Simulations of Expansion of a Plasma Cloud 

C B Chen and N Singh (Department of Electrical and Computer 
Engineering, University of Alabama in Huntsville, Huntsville, 
AL 35899, 205-895-6678) 

The dynamics and expansion of plasma clouds created by 
chemical releases in space is a complex problem. Some kinetic 
aspects of the dynamics are examined by means of one- and 
two-dimensional numerical simulations using Vlasov and PIC 
codes. The role of electron dynamics in the expansion is 
examined through the evolution of the electron velocity 
distribution function inside the cloud The two— dimensional 
PIC code is used to predict the potential structure when the size 
of the cloud is varied The evolution of the perpendicular 
electric fields near the boundaries of the cloud and their 
penetration into the ambient plasma affect the composite 
potential structure, which, in turn, affects the field-aligned 
plasma expansion Conditions under which field-aligned plasma 
expansion is the main feature of the dynamics of the cloud will 
be presented. 


SA51A-10 11201 

Generation of Equatorial Spread-F by Cosmoids - 
Comparison of Theory and Recent Experiments 

M Min (NASA/Goddard Space Flight Center, Code 916, Gretnbelt, 
MD 20771:301-286-5475) 

A theory for die formation of equatorial spread- F from the interaction 
in the F region of hypervelodry panicle streams formed from 
dispersed long-period meieoroids, coemoids, was presented several 
yens ago, but evidence for a cosmoid population was assumed and 
the upward motion of spread-F observed with incoherent scatter 
techniques could not be explained. The theory, however, was able to 
explain the formation of the electron bite-outs or bubbles, the pre 
midnight occurrence, the seasonal variations and latitude dependence 
offlpread-F. The theory explained the ion composition as the 
presence of Fe*, the fine structures observed with radar and the 
sdmUlaoons at GHz frequencies. Recent analysis of Pioneer 10 and 
1 1 measurements has supplied evidence for the existence of cosmoids 
In interplanetary space and their dispersion characteristics. A radar 
meteor signature of dispersed cosmoids has at last been detected at 2 
MHz to 140 km altitude. Chemical release experiments were used to 
demonstrate the formation of bite-outs by electron attachment and to 
detect Joule heating in the spread-F region. Shaped-charge chemical 
releases produced measurable ionization, demonstrating the ease of 
ion pair formation from critical velocity effects, thus explaining the 
upward motions measured with incoherent scatter radar. 

Furthermore, the discovery of sudden neutral sodium layers and 
sporadic Fe layers related to sporadic E represent another signature of 
the cosmoid population in the E region. An overview of the theory 
and experimental evidence will be presented. 


SA51A-11 113411 

Electric Field Formation at Solar Cell Edge Regions 

H Thiemann (Thiemann k Noack, Langenhardstr. 9, D-7800 
Freiburg. West Germany) 

R W Schunk (Center for Atmospheric and Space Sciences, Utah 
State University, Logan, Utah 84322-4405} 

Ground tests of negatively biased solar array test modules in 
a simulated LEO-plasma environment identified the solar cell 
edge region as a critical location. Arc discharges cause car- 
bonization spots along the cell peripheries at the K&pton sub- 
strate. Long-term exposure of solar arrays under such condi- 
tions will lead to short circuits and, eventually, a complete aolar 
array power low. For a quantitative assessment of this problem, 
we performed appropriate self-consistent computer experiments, 
simulating the solar cell-substrate- interface. The interaction is 
controlled by the electrical and the secondary emission proper- 
ties of the target. We will present results on the formation of 1 
strong electric fields, their dynamic features, strengths, direc- 
tions and locations. 


SA51A-12 1148(1 

The Dynamics of Finite-Sized, Low-fi, Plasma Clouds 
Injected at an Angle to a Background Magnetic Field 

T Neubgt and R H Miller (Both at: Space Physics Research 
Laboratory, University of Michigan, Ann Arbor, MI 48109- 
2143) 

O Buneman (Space, Telecommunications, and Radio science 
Laboratory, Stanford University, Stanford, CA 94305-4055) 

K-I Nishikawa (Department of Physics and Astronomy, 

University of Iowa, Iowa City, Iowa 52242) 

R Roussel -Dupre (Space Science and Technology, Los Alamos 
National Laboratory, Los Alamos, NM 87545) 

The dynamics of finite-sized, low-fi, plasma clouds moving at 
angles to a background magnetic field have been studied by means 
of a new three-dimensional electromagnetic panicle simulation 
code. The onset of electric and magnetic fields, the formation of 
charge layers, and the "braking" of the plasma cloud are studied for 
a range of ambient plasma and beam parameters. The electro- 
magnetic radiation for a finite- sized cloud is compared with earlier 
two-dimensional analytical and numerical results and to 
observations from sounding rockets, space shuttles and satellites in 
the Earth's ionosphere and magnetosphere 


SA52A CA: 313 Fri 1330h 

Chemical Release Experiments in Space II 

(joint with SM ) 

Presiding: P A Bernhardt, Naval 
Research Lab 

SA52A-1 1330h INVITED 

The NICARE Rxper^wMts -^^ pverviev 

Paul A Bernhardt (Space Plasma Branch, Plasna Physics 
Division, NavaS^ Research Laboratory, Washington, DC 
20375; 202-767-0196) 

The releases of electron attachment materials into the 
■id-latitude and equatorial ionospheric have been 
shown to be effective generators of (1) large scale 
plasma depletions, (2) enhanced plasma vaves, and (3) 
bright airglov clouds. During the Negative Ion, 
Cation Release Experiments (NICARE-1 and -2), nickel 
carbonyl, trif luromethyl bromide, and sulfur 
hexafluoride vere released between 300 and 380 km 
altitude into the dusk ionosphere. The experiments 

vere conducted from rockets launched from the Wallops 
Flight Facility, Virginia and from the Kvajalein 
Missile Range, Marshall Islands. For comparison, 

NICAFE-2 also contained a chemically inert, atomic 
release to estimate the effect of "snov-plov" density 
changes. The electron attachment molecules produced 
localized density reductions with factors ranging 
between 10 2 and 10 A . The density changes were 
measured with In situ Langmuir probes, radio beacon 
signals, and Incoherent scatter radar signals. Plasma 
waves from DC to over several MHz were stimulated by 
the releases. Some of the new wave modes can be 
directly traced to the presence of negative ions in 
the plasma. Enhanced optical emissions vere the 
result of chemiluminescent reactions between the 
neutral and plasma components of the releases and 
background environment. 

Research supported by ONR and NASA. 


SA52A-2 135® 

Simulations of High (Orbital) Velocity Releases of 
Electron Attachment Materials in the Ionosphere* 

W A Scales and P A Bernhardt (Space Plasma Branch, 
Plasma Physics Division, Naval Research Laboratory, 
Washington DC 20375; 202-767-3630) 

Two-dimensional numerical simulations of high (orbital) ve 
locity releases of electron attachment material clouds across 
the geomagnetic field in the F region are presented. The 
electrostatic simulation model includes electron attachment 
and mutual neutralization chemistry, self consistent electric 
fields, three species transport and ion inertia. The dvnam 
ics of the plasma irregularities produced are studied at early 
and at late times after release At early times, the plasma 
enhancement produced shows steeping of the backside and 
generation of vortices due to a Kelvin-Helmholtz instability 
At late times when the cloud has slowed due to collisions, 
the plasma depression produced shows structuring due to the 
E x B interchange instability if a background neutral wind 
is present. The effects of the variation of the ion-neutral col- 


lision frequency with altitude on the irregularity structuring 
will be discussed in detail 


•Research sponsored by ONT 


SA52A-3 14951) 

Electric Pield and Langmuir Probe Measurements in 
Chemically Modified Ionospheres* 

C L Si ef ring , P Rodriguez, D G Haas, M H Baumback, 

D N Walker , P A Bernhardt, (Space Plasma Branch, 
Plasma Physics Division, Naval Research 
Laboratory, Washington, DC 20375-5000) 

A sophisticated set of in-situ plasma diagnostics 
have been developed to study the perturbation 
effects of chemical releases on the ionosphere. 
Electric field, plasma vave and Langmuir probe data 
taken during the NICARE 16 2 experiments will be 
discussed with an emphasis on newly developed 
measurements. The Langmuir Probe measurements 
characterize the electron density depletions or 
enhancements and electric field measurements provide 
information on cloud polarization. Electrostatic 

plasma oscillations are almost always found in the 
highly disturbed regions near a chemical releases. 
Our plasma wave receiver can provide information on 
the spectral, temporal and k-space characteristics 
of the plasma vaves over a broad frequency range (DC 
to 10 MHz). Ve believe these measurements will 
answer many critical questions relating to 

effectiveness of various chemicals as ionospheric 
perturbation sources. 

♦Research supported by ONR and NASA. 


SA52A-4 1420b 

Negative Ion Plasma Vaves Observed in the NICARE 
Chemical Release Experiments* 

P Rodriguez , C L Siefring. D C Haas, P A Bernhardt, 
P V Nuygen, M M Baumback, D N Valker, (Space 
Plasma Branch, Plasma Physics Division, Naval 
Research Laboratory, Washington, DC 20375-5000) 

The release of several electron-attachment chemicals 
(SF6, Ni(C0)4, and CF3Br) in the ionosphere by the 
NICARE 1 and 2 sounding rockets have produced local 
plasmas in which the negative charge is carried by a 
heavy ion and in which there are none or very few 
electrons present. Such a plasma corresponds to a 
new environment for space experiments and ve have 
observed several varieties of low frequency 

electrostatic vaves which appear to be "negative ion 
plasma waves." Ve report the observations of waves 
based on electric field and plasma density 
measurements taken with the NRL electric field 
sensors and current probes. Electric field and 
plasma density fluctuations within the electron- 
depleted ionospheric "hole" are discussed, with an 
illustration of the evolution of the wave spectrum 
as the diagnostic vehicle leaves the region of 
electron depletion and re-enters the normal 

background ionosphere. The possible effects of wave 
normal orientations on the measured spectrum are 
discussed. It is believed that these experiments 
have produced the first collisionless space plasma 
in which negative ion plasma waves have been 
detected. „ 

♦Research supported by ONR and NASA. 


SA52A-5 14351) 

Far Ultraviolet Spectrometer Observations for the 
NICARE Rocket Experiment 

A B Christensen and J B Pranke (Both at: Space Sciences 

Laboratory, The Aerospace Corporation, PO Box 92957, Los 
Angeles, CA 90009) 

P A Bernhardt (Naval Research Laboratory, 4555 Overlook Avenue, 
5W, Washington, DC 20375-5000) 

A 1/8 m Eben-Fasuc spectrometer which scanned the far ultraviolet 
wavelength range from 125 to 149 nm at 1 .0 nm spectral resolution 
was included in the NICARE rocket payload launched 23 October 
1989 from Wallops Island, VA The instrument viewed 
perpendicular to the rocket spin axis and recorded emissions from OI 
(130.4), 01 (135.6) and a feature near 148.0 nm, possibly ft carbon 
line, that were enhanced during some of the chemical releases. The 
absolute intensity of the emissions and comparison with expected 
atrglow emission rates will be presented 


SA52A-6 14500 

HF Observations of Ionospheric Depletions 
Produced by Chemical Releases 

T. J. Fitzgerald. P. A. Argo, and R, C. 

Carlos (Los Alamos National Laboratory', 

SST-7. MS D466, P.O. Box 1663, Los Alamos, 

New Mexico 87545; (505) 667-1542) 

We have made observations of the ionospheric depletions 
produced by the NICARE l and 11 chemical releases using 
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Auroral Modeling of the 3371 A Emission Rate: 
Dependence on Characteristic Electron Energy 

P. G. Richards 

Computer Science Department and Center for Space Plasma and Aeronomic Research 
University of Alabama in Huntsvtlle 


D. G. Torr 


Phystcs Department and Center for Space Plasma and Aeronomic Research 
University of Alabama in Huntsville 


We have developed an efficient two-stream auroral electron model to study the deposition of 
auroral energy and the dependence of auroral emission rates on characteristic energy. This model 
incorporates the concept of average energy loss to reduce the computation time. Our simple two* 
stream model produces integrated emission rates that are in excellent agreement with the much 
more complex multistream model of Strickland et ol. (1983) but disagrees with a recent study 
by Rees and Lummerzheim (1989) that indicates that the N 2 second positive emission rate is a 
strongly decreasing function of the characteristic energy of the precipitating flux. Our calculations 
reveal that a 10 keV electron will undergo approximately 160 ionizing collisions with on average 
energy loss per collision of 62 eV before thermalizing. The secondary electrons are created with 
an average energy of 42 eV. When all processes including the backscattered escape fluxes are 
taken into account, the average energy loss per electron-ion pair is 35 eV in good agreement with 
laboratory results. 


1. Introduction 

There is currently renewed interest in the use of au- 
roral optical emission rates to deduce the characteris- 
tics of the precipitating particle fluxes, and ultimately, 
the global auroral energy input to the Earth’s upper 
atmosphere. Images from the Dynamics Explorer satel- 
lite have been used by Rees et al. [1988] to calculate 
the energetic electron flux and its characteristic energy. 
Imaging instruments planned for the ISTP mission will 
monitor key UV emissions on a global scale for the ex- 
press purpose of determining the global energy input. 

Early work in determining auroral particle character- 
istics from emissions concentrated on the use of the ra- 
tios of atomic oxygen emission rates (6300 A, 5577 A) to 
molecular nitrogen ion emission rates (3914 A, 4278 A) 
to deduce the incident auroral spectrum [Rees and 
Luckey, 1974; Vallance Jones , 1975; Shepherd et a/., 
1980; Strickland et al. , 1983]. The higher energy auroral 
electrons penetrate deeper into the thermosphere where 
the relative proportion of atomic oxygen is smaller. 
Thus the ratio of atomic to molecular emission rates de- 
creases with increasing electron energy. Unfortunately, 
chemical processes play an important role in the atomic 
oxygen emissions and it is difficult to separate the ef- 
fects caused by the characteristics of the auroral energy 
flux from the effects caused by changes in the atmo- 
spheric composition. Therefore, it would be useful to 
find an emission rate ratio that is sensitive to the au- 
roral characteristics but which is not complicated by 
chemical factors. 
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Recently, Rees and Lummerzheim [1989] suggested 
that the ratio of the second positive to first negative 
emission rates could be used to determine the charac- 
teristic energy of the auroral electron flux. Using an 
auroral electron model developed by Lummerzheim et 
al. [1989], Rees and Lummerzheim [1989] found that 
the N 2 second positive (3371 A ) emission rate de- 
creases substantially with increasing characteristic en- 
ergy of the auroral electrons while the Nt emission 
rates are almost constant. This ratio would be an at- 
tractive alternative to those used previously because it 
would be independent of atmospheric composition and 
both emissions are prompt, thus eliminating chemical 
effects. Unfortunately, the calculations of Rees and 
Lummerzheim [1989] are in conflict with the earlier cal- 
culations by Daniell and Strickland [1986] who found 
that the 3371 A emission rate was nearly independent 
of the characteristic energy. 

The experimental evidence also seems to be in con- 
flict. Rees and Lummerzheim [1989] present data from 
high flying aircraft that support their theoretical calcu- 
lations. On the other hand, Solomon [1989] presented 
data from the visible airglow instrument on the At- 
mosphere Explorer C satellite showing that the ratio 
of the N 2 3371 A to 4278 A emission rates has 
only a small dependence on the characteristic energy, 
which can be accounted for by contamination of the 
3371 A second positive emission by the Vegard-Kaplan 
(0-9) band. The VAE data support the earlier calcula- 
tions of Daniell and Strickland [1986] and Strickland et 
al. [1983]. Solomon was able to reproduce the observed 
ratios using his own two-stream auroral electron depo- 
sition code. We note that the experimental data pre- 
sented by Solomon [1989] for the ratio of N 2 3371 A to 
N i 4278 A is in excellent agreement with the ratio of 
N 2 3371 A to N+ 3914 A that was measured on a 1974 
rocket flight by Sharp et al. [1979]. 



PRECEDING PAGE BLANK NOT FILMED 


i 0.338 Richards and Torr: Auroral 3371 A Emission Rate 


The source of the discrepancy between the models is 
difficult to understand. Both Rees and L ummerzheim 
[1989] and Strickland et ai [1983] use relatively com- 
plete multi-stream models. The model used by Solomon 
[1989] is a simpler two-stream model and his calcula- 
tions are in agreement with the calculations of Strick- 
land et ai [1983] and Dameli id Str-ckland [1986] 
which show little variation in the 3371 A to 
3914 A emission rate ratio as a function of characteristic 
energy. In this paper we present calculations of a num- 
ber of important emission rates obtained from a two- 
stream auroral electron model that we have developed. 
The calculated emission rates for the 3371 A emission 
rate and a number of other emission rates are in accord 
with the earlier calculations of Strickland et al. [1983] 
and Danxell and Strickland [1986] for the dependence 
on characteristic energy of the precipitating flux. How- 
ever, mainly due to the use of a revised cross section, 
our atomic oxygen 1356 A emission rate is a factor 
of 2.5 lower than that calculated by Strickland et ai 
[1983]. 

2. Model 

2.1. General Principles 

The model that we have developed is based on the 
two-stream photoelectron flux model of Nagy and Banks 
[1970] that was subsequently extended to 500 eV and 
combined with a continuous energy loss model to calcu- 
late auroral electron fluxes by Banks et ai [1974]. Al- 
though both our model and the model of Solomon et ai 
[1988] have origins in the Nagy and Banks [1970] two- 
stream model, they have evolved substantially along en- 
tirely separate paths. Our model owes much to our ear- 
lier work with the ionospheric photoelectron flux 
[Richards and Torr , 1984; 1985a]. 

By incorporating a variable energy grid developed by 
Swartz [1985] and a variable altitude grid we have been 
able to extend the two-stream model up to energies 
greater than 20 keV. This has eliminated problems en- 
countered by Banks et ai [1974] in matching the con- 
tinuous slowing down approach that they used above 
500 eV with the two-stream approach they used below 
500 eV. The continuously variable altitude grid allows 
altitude steps of less than a kilometer below 90 km up 
to 50 km at 500 kin with a manageable number of grid 
points. 

Further economy in computer time is achieved by in- 
troducing the concept of an average energy loss for an 
excitation or an ionization event. This concept allows 
the use of only the total excitation and ionization cross 
sections instead of treating each partial excitation or 
ionization process separately. That is, the excitation 
(or ionization) is treated as arising from a single aver- 
age state for each thermospheric species. The average 
energy loss depends on the species and also on the en- 
ergy of the primary electron. The calculation of the 
emission rates then becomes a two stage process with 
just the total cross sections being needed to calculate 
the electron flux as a function of energy and altitude 
in the first stage. In the second stage, the partial cross 
sections are folded with the electron fluxes to produce 
the excitation rates. 


For each of the 3 main neutral atmospheric compo- 
nents (O, O 2 , N 2 ) the total cross sections are made up 
of partial cross sections from the numerous electronic 
states of each species each corresponding to a different 
energy state of the atom or molecule. In addition, each 
electronic state of a molecule can be created in any one 
of a large number of vibrational energy states. A com- 
plete evaluation of the electron flux would require the 
separate accounting of all these energy losses which we 
have replaced with a single, energy dependent average 
energy loss for the excitation and ionization of each of 
the three main thermospheric species O, 0 2 , and N 2 . 

The treatment of excitation processes is relatively 
straightforward; the average energy loss is specified and 
the electrons deposited in the correct lower energy bin. 
However, ionizing collisions are more complex because 
of the production of secondary electrons which may be 
produced with energies ranging from 0 up to E v - I{ 
where E v is the energy of the primary electron and U is 
the ionization potential of the state i being produced. 
We follow the approach of Banks et ai [1974] and des- 
ignate the higher energy electron as the degraded pri- 
mary and the lower energy electron as the secondary. 
This means that the maximum secondary energy is then 
( E p — /i)/2 as is the minimum energy of the degraded 
primary. We treat ionization as arising from a single 
state with an average ionization potential / and we use 
the measured secondary electron distributions of Opal 
et ai [1971] to determine the average energy E f of the 
secondary electrons produced by an electron of energy 
E v (we note that the measured secondary electron dis- 
tributions are in fact a sum of the contributions from 
all the ionization states). The average energy of the 
degraded primaries is then E v — I — E A separate 
ionization potential could be used for each species but 
a further improvement in computational efficiency can 
be made by observing that the ionization potentials and 
secondary electron distributions are similar enough that 
a single ionization potential and secondary electron dis- 
tribution for all 3 major species will suffice. Since N 2 is 
the most important species in aurora, we adopt the N 2 
ionization characteristics. At each electron energy, the 
total number of secondary electrons from O, 0 2 , and N 2 
is calculated and then they are distributed according to 
the measured secondary electron distributions of Opal 
et ai [1971]. The model has been found to conserve 
energy to better than 5% for characteristic energies in 
the range .1 to 20 keV. 

2.2. Average Ionization Potentials 

The calculation of the average ionization potential 
follows from the knowledge of the ionization potentials 
and the partial cross sections for the various ionization 
states of each of the thermospheric species. For exam- 
ple, N 2 ionization results in the formation of the X, A, 
and B states of Nt and also N + which arises from sev- 
eral higher lying states that dissociate. The ionization 
potentials for the X, A, and B states are 15.6, 16.8, and 
18.8 eV respectively while the bulk of the N + probably 
arises from a state with a threshold near 37 eV [ Erd - 
man and Zip } [ 1986]. The actual energy lost by the 
primary electrons may be greater than threshold due 
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to the ions being vibrationally excited and because, in 
the case of dissociative ionization, the atomic fragments 
are observed to carry substantial kinetic energy [Lil- 
iana and St or dale , 1975]. To complete the calculation 
of t lie average ionization energy for Ni we also need 
the ratios of the partk*. cross sections. We examine 
the high energy case first where the proportions of the 
various partial cross sections are practically constant 
and the calculation is straightforward. Above approxi- 
mately 100 eV, dissociation accounts for —20% [Rapp et 
al. , 1965] and the B state —10% [Borst ar\d Zipf , 1970] 
of the total cross section. According to Cartwright et 
al. [1975], the X and A state cross sections are ap- 
proximately equal (that is, —35%). Using these per- 
centages for the partial cross sections and ionization 
potentials given above, the average ionization potential 
is —20 eV for high energy electron impact ionization of 
No. We now examine the low energy average energy loss 
per collision. The average ionization potential begins 
at 15.6 eV at threshold when the X state is produced 
but increases to —16 eV above 17 eV when the A stat'e 
threshold is reached. There is only a marginal increase 
when the B state threshold is reached at 18.8 eV because 
it accounts for less than 10% of the total cross section 
at these energies and the energy loss of 18.8 eV is only 
marginally larger than 16.8 eV. After about 30 eV dis- 
sociative ionization becomes important. The ionization 
potential then rises steadily to approximately 19 eV at 
50 eV before leveling off toward 20 eV at high energies. 

A similar calculation yields an average ionization po- 
tential of approximately 18 eV for both O and O 2 , at 
high energies. Because N 2 is the dominant constituent 
and all three ionization potentials are similar, we use 
he average N 2 ionization potential for all three species. 
This simplifies the calculation and reduces the compu- 
tation time without introducing significant errors. The 
adopted ionization and excitation energy losses as a 
function of primary electron energy are shown in Fig- 
ure 1. This figure also shows the average energy of the 
secondary electrons as a function of the primary en- 
ergy. Below 25 eV the ionization potential is assumed 
to be 16 eV while above 25 eV it is represented by 
10(1 + (1 - 15 /E) 112 ) eV. 


2.3. Average Secondary Electron Energy 

To obtain the secondary electron distribution we 
adopt the Banks et al. [1974] parameterization of the 
Opal et al. [1971] measured distributions. The proba- 
bility of a secondary electron of energy E, is given by 


P t (E.) 


A 

1 + {EJEY 


(1) 


where A = [(£* arctan (£ m / E)~ 1 } is a normalization fac- 
tor that ensures a total probability of unity when in- 
tegrated over all secondary energies from 0 to E in = 
( E v — I )/ 2 and E=14 is an empirical normalization fac- 
tor. The product of E , and Equation 1 integrated over 
energy yields the average secondary electron energy for 
a primary energy E v as 


( 2 ) 



Fig. 1. Average energy losses per collision and average secondary 
electron energies (E„«) as a function of primary electron energy. 
The average ionization potential is labelled I and the total en- 
ergy loss per ionization is labelled I+E U t, The average excitation 
potentials are indicated by an asterisk and the Oj excitation po- 
tential is set equal to that of Nj . Note that below 5 eV, the N 2 
excitation potential is set at 1 eV. 


The dependence of the average secondary energy on 
the primary energy is shown in Figure 1 along with 
the average total energy loss of the primary (I + £**)• 
The average energy of the secondaries increases steadily 
from 0 near threshold to 52 eV for 10 keV primary elec- 
trons. The total energy loss per ionizing collision for 
10 keV electrons is 72 eV, when the 20 eV ionization 
potential is included. The average energy of the secon- 
daries is approximately equal to the ionization energy 
for 200 eV primary electrons. 

A quantity of interest, in relation to energy degra- 
dation of high energy electrons, is the average energy 
required to produce each electron-ion pair. The aver- 
age energy required to produce eich electron-ion pair 
is a quantity that is independent of electron energy 
and is also remarkably independent of the species be- 
ing ionized. The experimental value for the energy lost 
per electron-ion pair for high energy electron is 35 eV 
[Valentine and Curran , 1958]. This energy loss per 
electron-ion pair was used in early auroral electron de- 
position codes [Rees, 1963; Rees et al. 1969; Rees and 
Jones, 1973]. 

An approximate value for the energy loss per electron- 
ion pair can be deduced by using the average energy 
losses depicted in Figure 1, assuming that for electron 
energies above approximately 100 eV the energy lost 
to excitation collisions is small and can be neglected. 
With this assumption, it requires about 160 ionizing 
collisions to thernialize a single 10 keV electron. Thus, 
on the average, 62 eV is lost in the creation of each of 
the 160 electron-ion pairs. Since the ionization energy 
is 20 eV for electron energies above 100 eV, the average 
energy of the secondaries is —42 eV. This means that 


E„= 0.5 A E 2 hx[l + (EJE) 2 ) 
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the 160 secondary electrons are sufficiently energetic to 
create one more electron-ion pair each. If all secon- 
daries created an additional electron-ion pair, the total 
number of pairs would be 320 and the average energy 
per electron-ion pair would decrease from 62 to 31 eV. 
In reality this does not happen because, below 100 eV, 
excitation processes begin to compete effectively with 
the ionization processes for the available electron en- 
ergy and the number of additional electron-ion pairs 
produced by the secondary electrons would be less than 
160. In fact, for a 42 eV electron the total ionization 
and excitation cross sections are approximately equal 
and only half the secondaries could be expected to pro- 
duce an additional electron-ion pair. This agrees with 
our previous calculations that show that a third rather 
than a half, of the total ionization is created by de- 
graded primaries and secondaries with energies below 
100 eV [Richards and Tbrr, 1985b]. Thus the original 
10 keV electron would ultimately produce about 240 
electron ion pairs and yield an average energy loss per 
electron-ion pair of 42 eV. This is only an estimate of 
the energy lost per electron-ion pair and a more detailed 
calculation including transport is required to determine 
the actual value. It was pointed out by Banks et ai 
[1971] that escaping backscattered electrons will be lost 
to the system and act to increase the energy loss per 
electron-ion pair. 

We have summed the total ion production rate in our 
full auroral calculation, and we obtain an average en- 
ergy loss per electron-ion pair of 35 eV, which is smaller 
than our estimate but in agreement with the laboratory 
measured value. A slightly higher value of 37 eV was 
obtained by Fox and Victor [1988] using their discrete 
local energy loss method. The reason that the energy 
lost per electron-ion pair is not a strong function of 
electron energy has to do with the relationship between 
the average secondary electron energy and the primary 
electron energy. Electrons with higher initial energies 
suffer a greater energy loss per collision as they degrade 
but they produce higher energy secondaries which are 
more likely to generate secondary ionization. For ex- 
ample, a 1 keV electron will undergo only 22 ionizing 
collisions with an average energy loss of 45 eV before 
it thermalizes. Thus, the secondary electrons have an 
average energy of only 25 eV compared to the 42 eV for 
the 10 keV electrons and are much less likely to produce 
additional ions. 

2.4. Cross Sections 

In a number of previous studies of the ionospheric 
photoelectron flux we have chosen measured total exci- 
tation cross sections that produce good agreement be- 
tween theory and the photoelectron spectrometer mea- 
surements from the AE-E satellite [Lee et a/. t 1980]. At 
low altitudes where N 2 is the dominant species, the to- 
tal cross section obtained from electron mobility studies 
by Pitchford and Phelps [1982] was found to be com- 
patible with the PES measurements. At high altitudes 
where atomic oxygen is the dominant species, the emis- 
sion cross sections measured by Zipf and co-workers pro- 
duced good agreement between theory and PES mea- 
surements. These total cross sections for energies be- 


low 100 eV have been published by Richards and Torr 
[1988]. Basically the same cross sections have been used 
in this study but they have been reparameterized to ex- 
tend them to higher energies. 

Above 100 eV, the excitation cross sections decay 
rapidly with increasing energy and are much less im- 
portant than the ionization cross sections both because 
they are smaller and because the energy loss per colli- 
sion is smaller. The total ionization cross sections are 
better established than the total excitation cross sec- 
tions although there are some differences [Kieffer and 
Dunn , 1966]. We have adopted the N 2 and 0 2 total 
ionization cross sections of Rapp and Englander- Golden 
[1965] which have also been used by most other mod- 
elers. The total ionization cross section for O is from 
Brook et al. [1978], The elastic cross sections are very 
important because of their role in inhibiting transport. 
We have used the elastic cross sections of Solomon et 
al. [1988] and also their elastic backscatter coefficients. 
The total cross sections used in our auroral model are 
shown in Figure 2. Our N 2 total excitation cross section 
is comparable to that of Solomon et al. [1988] below 
25 eV but is smaller at higher energies. The differences 
at high energies have little effect on the calculated fluxes 
because the excitation cross section is smaller than the 
ionization cross section. However, differences in cross 
sections below 30 eV produce comparable differences in 
fluxes. The N 2 total excitation cross section of Strick- 
land et al. [1983] is almost a factor of two larger than 
ours at all energies and their fluxes would be a fac- 
tor of two lower below 30 eV, at least below 200 km 
where N 2 is the dominant species. The cross sections 
of Solomon et al. [1988] and Strickland et al. [1983] 
were obtained by summing the partial cross sections 
and there is the possibility of double counting some 
cross sections; for example, those that lead to dissocia- 
tion. Moreover, Strickland et al. [1983] included large 



Fig. 2. Totnl elastic, excitation, and ionization cross sections em- 
ployed in the model. The ionization cross sections are indicated 
by a plus while the excitation cross sections are indicated by an 
asterisk. 
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Rydberg cross sections from Green and Stolarski [1972] 
which have been revised sharply downward by Porter 
et al. [1976]. There is now good agreement between 
the total cross section of Pitchford and Phelps [1982] 
and the sum of the partial excitation cross sections of 
Cartwright et al. [1977a,b] as revised by Trajmar et al. 
[1983], below 20 eV. Above 20 eV, the ionization cross 
section becomes an increasingly important component 
of the total inelastic cross section and it is not easy to 
compare the two cross sections. 

The total atomic oxygen excitation cross section em- 
ployed by Solomon et al. [1988] is a factor of 2 larger 
than ours above 15 eV and will produce a similar dif- 
ference in flux above 250 km where O is the dominant 
species but the atomic oxygen cross section has little 
effect on the integrated emission rates. Our atomic oxy- 
gen excitation cross section was obtained by summing 
the measured emission cross sections for 1304, 1356, 
and 1027 A [. Zipf and Erdman, 1985], the 989 cross 
section from Gulcicek and Doering [1988], and the the- 
oretical 1 D and l S from Henry et al. [1969]. Implicit in 
this procedure is the assumption that the higher lying 
triplet and quintet states are included in the 1304 and 
1356 emission cross sections via cascade. We have left 
out some theoretical Rydberg cross sections proposed 
by Jackman et al. [1977] and some minor states that 
radiate directly to the ground state but for which there 
is no experimental data. Thus, our cross section must 
be regarded as a lower limit. 

Figure 3 shows the excitation cross sections for the 
second positive (C°T, t ) and Lyman-Birge-Hopfield 
(a 1 *-,,) systems of N 2 . Also shown is the cross section 
used for calculation of the 0( 5 5) 1356 A emission rate. 
We obtained this cross section by reducing the measured 
cross section of Stone and Zipf [1974] by the factor 3.1 
which is the same factor that the 1304 A emission cross 
section of Stone and Zipf [1974] was reduced by Zipf and 
Erdman [1985]. The (C 3 7t„) cross section was obtained 



Fig. 3. Cross sections for the three excited states giving rise to 
the emissions studied in this paper. 


by multiplying the 3371 A cross section of Imami and 
Borst [1974] by 4 and the (a l x,j) cross section is from 
Ajello and Shemansky [1985]. 

We have examined the sensitivity of the emission rate 
ratios to cross sections and this will be discussed later. 
In all these calculations we have used the 1 erg cm~ 2 s _1 
Gaussian incident flux distribution, and the neutral at- 
mosphere employed by Strickland et al. [1983]. 


3. Results 


3.1. Comparison With Previous Work 

We have calculated the N 2 3371 A , Nj 3914 A , Nj 
4278 A , O 1356 A , and several N 2 LBH band emission 
rates as a function of energy and these are shown in 
Figure 4. This figure shows that both the 3371 A and 
3914 A emission rates are independent of the char- 
acteristic energy of the precipitating flux for energies 
above 2 keV in agreement with the results of Strickland 
et al. [1983] and Daniell and Strickland [1986]. Not 
only is the shape in good agreement but, except for the 
1356 A emission rates, the magnitudes are also in good 
agreement. Although the shape of the 1356 A curve is 
in good agreement with that of Strickland et ai [1983], 
the magnitude is a factor of 2.5 lower owing to the use 
of the revised cross section of Zipf and Erdman [1985]. 



ENERGY (keV) 


Fig. 4. Circulated emission rates as a function of the characteris- 
tic energy for a Gaussian energy distribution with a total incident 
energy flux of 1 erg cm” 2 s” 1 . When differences in cross sections 
are taken into account, there is exceUent agreement with the cal- 
culations in Figure 8 of Strickland et al. (1983) and Figure 11 of 
Daniell and Strickland (1986). 
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We have included O 2 Schuman-Runge absorption which 
affects both the 1356 A and the 1200-1600 A LBH 
bands when the characteristic energies are high and the 
electrons penetrate to lower altitudes [Strickland et al . , 
1983]. The ratios of the LBH bands are taken from 
Ajello and Shemansky [1985]. 

Our calculated second positive to first negative emis- 
sion rate ratios are within 20% of the measured val- 
ues for the Nt emissions. At 10 keV, our 3371 A to 
3914 A ratio is 0.3 compared to 0.25 from Sharp et al. 
[1979] and our 3371 A to 4278 A ratio is 0.98 compared 
to 0.8 from Solomon [1989]. Solomon obtained better 
agreement between his model ratios and the measured 
ratios but there is sufficient uncertainty in the input 
parameters to account for the differences. 

3.2. Sensitivity of Ratios 


1 



UJ 



We have performed some parameter studies to char- 
acterize the sensitivity of the ratio of the second positive 
to first negative integrated emission rates to possible 
errors m the model inputs. Obviously, a reduction of 
20% in the 3371 A emission rate cross section would 
bring the calculated and measured values into excellent 
agreement but a 30% increase in the No total excitation 
cross section has a similar effect by decreasing the low 
energy electron flux which is responsible for most of the 
3371 A emission. Likewise, a 30% increase in the N 2 
total ionization cross section above 100 eV reduces the 
ratio from 0.3 to 0.25 by increasing the 3914 A pro- 
duct n rate. The integrated ratio is not sensitive to 
changes of up to a factor of two in most other param- 
eters including: the atomic and molecular oxygen in- 
elastic cross sections, the O, O 2 , and No elastic cross 
sections and backscatter coefficients, and the relative 
concentrations of the species. We estimate a possible 
error of 10% in our computed average excitation and 
ionization potentials but this has negligible effect on 
our computed ratios. 


ENERGY (ktV) 

Fig. 5. The initial partitioning of the incident 1 erg cm -3 s*** 1 
energy flux between ionization, excitation, thermal electron heat- 
ing, and backscatter as a function of characteristic energy. The 
largest proportion of the energy (^-35%) goes in::ially into the 
ionization potential of the while (^20%) goes into N? excita- 
tion. Only (^16%) is backscattered out of the thermosphere. O 
is an important absorber of energy at the lowest energies while 
becomes increasingly important as the characteristic energy 
ncreascs. 


3.4. Electron Flux Spectra 

Downward moving fluxes at 120, 174, 223, and 326 
km are shown in Figure 6 for a 5 keV incident Gaussian 
flux with an energy flux of 1 erg cm“ 2 s“ 1 . The incident 
flux can be seen centered at 5 keV in Figure 6. At the 
two highest altitudes, there is very little degradation 
of this initial flux but the degradation is noticeable at 
174 km and pronounced at 120 km. Because there is so 


3.3. Energy Budget 

The incident electron energy flux is initially parti- 
tioned into a large number of excitation and ionization 
processes before it finally emerges as heat for the ther- 
mosphere or is radiated into space. Figure 5 shows the 
gross energy partition amongst ionizations, excitations, 
thermal iectron heating and backscattered escape flux 
as a function of characteristic energy. Nt ions capture 
the greatest share of the available energy (35%). Exci- 
tation of No is next with (20%) while only 15% is lost 
through he escaping backscattered flux. This escape 
flux is ch smaller than the 45% obtained by Banks 

et al. i *4], possibly as a result of the use of different 

cross sec: : 011 s and backscatter coefficients. Below 1 keV, 
' uization and excitation of atomic oxygen absorption 

e important energy sinks for the electron energy, but 
they become small for high energy incident fluxes. Ab- 
sorption in molecular oxygen shows the opposite trend, 
becoming more important with the deeper penetration 
of the higher energy fluxes. Thermal electrons capture 
a greater proportion of the available energy, the lower 
the characteristic en rgy. 



ENERGY («V) 

Fig. 6. Downward moving electron flux spectra at several alti- 
tudes for a 5 keV Gaussian incident flux. The incident energy 
flux is 1 erg cm~ 3 s~* . 
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little interaction with the thermosphere for high energy 
electrons at the high altitudes, there are few degraded 
primaries to fill in the region between 300 eV and 5 keV. 
However, at the lowest altitudes this intermediate en- 
ergy range is filled in. Comparison of the downward 
fluxes in Figure 6 with the upward fluxes in Figure 7 
reveals that, below 225 km, where transport is inhib- 
ited, the electron flux is isotropic for energies less than 
300 eV. At 326 km, the upward (escape) flux is a factor 
of 2 larger than the downward flux at low energies and 
orders of magnitude larger at intermediate energies. 

4. Conclusions 

We have developed an efficient two-stream auroral 
electron model that incorporates the concept of aver- 
age energy loss. This model produces integrated emis- 
sion rates that are in excellent agreement with the more 
sophisticated multi-stream model of Strickland et ai 
[1983] but is in disagreement with the model of Rees and 
Lummerzheim [1989] with regards to the energy depen- 
dence of the Nt 3371 A second positive emission rate. 
Our calculations give a value of 35 eV for the average 
energy lost per electron-ion pair produced independent 
of primary electron energy and we have explained this 
behavior in terms of the variation in the energy of the 
secondary electrons. We find that more than 30% of 
the initial energy flux is stored initially as ionization 
energy of Nj while about 20% goes into excited states 
of Nt while only 15% is backscattered out of the ther- 
mosphere. All other processes are minor except at low 
incident energies where 20% of the energy is stored in 
atomic oxygen ions. 
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Fig. 7. Upward moving electron flux spectra for the 5 keV Gaus- 
sian incident flux shown in Figure 6. 


References 

Ajello, J. M., and D. E. Shemansky, A reexamination 
of important N 2 cross sections by electron impact 
with application to the day glow: The Lyman- Birge- 
Hopfield band system and N I (119.99 nm), J. Geo - 
phys. Res., 90, (A10), 1985. 

Banks, P. M. and C. R. Chappell, and A. F. Nagy, 
A new model for the interaction of auroral electrons 
with the atmosphere: spectral degradation, backscat- 
ter, optical emission, and ionization, J. Geophys. Res., 
79, 1459, 1974. 

Borst, W. L., and E. C. Zipf, Cross section for electron- 
impact excitation of the (0,0) first negative band of 
Nj from threshold to 3 keV, Phys. Rev. A , 1 , 834, 
1970. 

Brook, E., M. F. A. Harrison, and A. C. H. Smith, 
Measurements of the electron impact ionisation cross 
sections of He, C, O and N atoms, J . Phys. B, 11 , 
3115, 1978. 

Cartwright, D. C., W. R. Pendleton, and L. D. Weaver, 
Auroral emission of the Nj Meinel bands, J. Geophys. 
Res., SO, 651, 1975. 

Cartwright, D. C., A. Chutjian, S. Trajmar, and W. 
Williams, Electron impact excitation of the electronic 
states of N 2 , I, Differential cross sections at incident 
energies from 10 to 50 eV, Phys. Rev . A, 16, 1013, 
1977a. 

Cartwright, D. C., S. Trajmar, A. Chutjian, and W. 
Williams, Electron impact excitation of the electronic 
states of N 2 , II, Integral cross sections at incident 
energies from 10 to 50 eV, Phys. Rev. A , 16, 1041, 
1977b. 

Daniell, R. E., and D. J. Strickland, Dependence of the 
auroral middle UV emissions on the incident electron 
spectrum and neutral atmosphere, J. Geophys. Res., 
91, 321, 1986. 

Erdman, P. W., and E. C. Zipf, Dissociative excitation 
of the NI( 5 S) state by electron impact on N 2 : Exci- 
tation function and quenching, J. Geophys. Res., 91 , 
1986. 

Fox, J. L., and G. A. Victor, Electron energy deposition 
in N 2 gas, Planet. Space Sci., 36, 326, 1988. 

Garstang, R. H., Multiplet intensities for lines 4 S — 2 D 
of SII, Oil and NI, Astrophys . J., 115, 506, 1952. 

Green, A. E. S. and R. S. Sfcolarski, Analytic models of 
electron impact excitation cross sections, J. Atmos. 
Terr. Phys., 3 4, 1703, 1972. 

Gulcicek, E. E., and J. P. Doering, Absolute differential 
and integral electron excitation cross sections of the 
atomic oxygen 3 P and 5 P states at 30 eV, J. Geophys. 
Res., 92{ A4), 3445-3448, 1988. 

Henry, R. J. W., P. G. Burke and A. L. Sinfailam, Scat- 
tering of electrons by C, N, O, N + , 0+ and 0 ++ , 
Phys. Rev., 178 218, 1969. 

Jackman, C. H., R. M. Garvey, and A. E. S. Green, 
Electron impact on atmospheric gases, I. Updated 
cross sections, J. Geophys. Res., 82, 5081-5090, 1977. 



10.344 


Richards and Torr: Auroral 3371 A Emission Rate 


Kieffer, L. J., and G. H. Dunn, Electron impact ion- 
ization cross section data for atoms, atomic ions and 
diatomic molecules: I. experimental data, Rev. Mod . 
Phys., 38 , 1, 1966. 

Lee, J. S., J. P. Doering, T. A. Potemra, and L. H. 
Brace, Measurements of the ambient photoelectron 
spectrum from Atmosphere Explorer, I, AE-E mea- 
surements below 300 km during solar minimum con- 
ditions, Planet. Space Set., 28, 947, 1980. 

Liliana, D, and J. A. D. Stockdale, Dissociative ioniza- 
tion of molecules by electron impact, II, Kinetic en- 
ergy and angular distributions of N+ and N ++ ions 
from N;, J . Chem. Pkys 63 , 3898-3906, 1975. 

Lununerzheim, D., M. H. Rees, and H. R. Anderson, 
Angular dependent transport of auroral electrons in 
the upper atmosphere, Planet . Space Sci., 37, 109, 
1989, 

Malunood, I, and Walter L. Borst, Electron excitation 
of the (0,0) second positive band of nitrogen from 
threshold to 1000 eV*, J. Chem. Phys ., 61 , 1115, 
1974. 

Nagy, A. F., and P. M. Banks, Photoelectron fluxes in 
the ionosphere, J. Geophys. Res., 75, 6260, 1970. 

Opal, C. B., W. K. Peterson, and E. C. Beaty, Mea- 
surements of secondary-electron spectra produced by 
electron impact ionization of a number of simple gases, 
/. Chem. Phys., 55, 4100, 1971. 

Pitchford, L. C., and A. V. Phelps, Comparative calcu- 
lations of electron swarm properties in at moder- 
ate E/N values, Pkys. Rev. A 25 , 540, 1982. 

Porter, H. S., C. H. Jackman, and A. E. S. Green, Effi- 
ciencies for production of atomic nitrogen and oxygen 
by relativisitic proton impact in air, J. Chem. Phys., 
65, 154-167, 1976. 

Rapp, D., and P. Englander-Golden, Total cross sec- 
tions for ionization and attachment in gases by elec- 
tron impact I, Positive ionization. J. Chem. Phys., 
43, 1464, 1965. 

Rees, M. H., Auroral ionization and excitation by inci- 
dent energetic electrons, Planet. Space Sci., 11 , 1209, 
1963. 

Rees, M. H. and R. A. Jones, Time dependent studies 
of the aurora, II, Spectroscopic morphology, Planet. 
Space Sci., 21, 1213, 1973. 

Rees, M. H., and D. Luckey, Auroral electron energy 
derived from the ratio of spectroscopic emissions, 1, 
Model computations, J. Geophys. Res., 79, 5181, 
1974. 

Rees, M. H., and D. Lummerzheim, Characteristics of 
auroral electron precipitation derived from optical spec- 
troscopy, J. Geophys. Res., 94, 6799, 1989. 

Rees, M. H., D. Lummerzheim, R. G. Roble, J. D. Win- 
liinghani, J. D. Craven, and L. A. Frank, Auroral En- 
ergy deposition rate, characteristic electron energy 
and ionospheric parameters derived from Dynamics 
Explorer 1 images, J. Geophys. Res., 93, 12841, 1988. 

Rees, M. H., A. I. Stewart, and J. C. G. Walker, Sec- 


ondary electrons in aurora, Planet. Space Sci., 17, 
1997, 1969. 

Richards, P. G., and D. G. Torr, An investigation of the 
consistency of the ionospheric measurements of the 
photoelectron flux and solar E"V flux, J. Geophys. 
Res., 89, 5625, 1984. 

Richards, P. G., and D. G. Torr, The altitude variation 
of the ionospheric photoelectron flux: a comparison 
of theory and measurement, J. Geophys. Res., 90, 
2877, 1985a. 

Richards, P. G. and D. G. Torr, On the production 
of N + by energetic electrons, J. Geophys. Res. 90, 
9917, 1985b. 

Richards, P. G., and D. G. Torr, Ratio of photoelec- 
tron to EUV ionization rates for aeronomic studies 
J. Geophys. Res., 93, 4060, 1988. 

Sharp, W. E., M. H. Rees, and A. I. Stewart, Coordi- 
nated rocket and satellite measurements of an auro- 
ral event, 2, The rocket observations and analysis, J. 
Geophys. Res., 84, 1977, 1979. 

Shepherd, G. G., J. D. Winningham, F. E. Bunn, and 
F. W. Thirkettle, An empirical determination of the 
production efficiency for auroral 6300- A emission by 
energetic electrons, J. Geophys . Res., 85, 715, 1980. 

Solomon, S. C., Auroral excitation of the N 2 2P(0,0) 
and VK(0,9) bands, J. Geophys. Res., 94, 17215, 
1989. 

Solomon, S. C., P. B. Hays, and V. J. Abreu, The au- 
roral 6300 A emission: Observations and modeling, 
J. Geophys. Res., 93( A9), 9867, 1988. 

Stone, E.J., and E.C. Zipf, Electron-impact excitation 
of the 3 S° and s S p states of atomic oxygen, J , Chem. 
Phys., 60, 4237, 1974. 

Strickland, D. J., J. R. Jasperse, and J. A. Whalen, De- 
pendence of auroral FUV emissions on the incident 
electron spectrum and neutral atmosphere, J. Geo- 
phys. Res., 88, 8051, 1983. 

Swartz, W. E., Optimization of energetic electron en- 
ergy degradation calculations, J. Geophys. Res., 90, 
6587, 1985. 

Trajmar, S., D. F. Register, and A. Chutjian, Electron 
scattering by molecules, II, Experimental methods 
and data, Phys. Rep., 97, 219, 1983. 

Vallance Jones, A., A model for the excitation of elec- 
tron -urora and some applications, Can . J. Phys., 
53, 2276, 1975. 

Zipf, E. C\, and P. W. Erdman, Electron impact exci- 
tation of atomic oxygen: Revised cross sections, J. 
Geopnys. Res., 90, 11087, 1985. 


P. G. Richards and D. G. Torr, University of Alabama, Re- 
search Institute C-10, Huntsville, AL 35899. 


(Recieved October 18, 1989; 
accepted November 28, 1989.) 



CR'G^AL PAGE IS 

OF POOR QUALITY 


JOURNAL OF GEOPHYSICAL RESEARCH. VOL. 95. NO. A6. PAGES 77:5-7733. JUNE I. 1990 


The Dependence of Modeled 01 1356 and Lyman Birge Hopfield 
Auroral Emissions on the Neutral Atmosphere 

G. A. Germany and M. R. Torr 

Space Science Laboratory. \'ASA Marshall Space h'livht Center Huntsville. Alabama 

P. G. Richards and D. G. Torr 

University of Alabama in Huntsville. Huntsville . Alabama 


Images of the enure auroral oval al carefully selected wavelengths contain information on the global energy 
influx due to energetic particles and some information on the characteristic energy of the precipitating particles 
In this paper we investigate the sensitivity of selected auroral emissions to changes m the neutral atmosphere. In 
particular, we examineW behavior of 01 1356 A and two Lyman Birge Hopfield (LBH) bands and their ratios 
to each other with changing atmospheric composition The two LBH bands are selected so that one lies in the 
region of strong CL absorption 1 1464 A I and one lies at a wavelength where CL absorption is effectively negli- 
gible f 1838 A). We find that for anticipated average uncertainties in the neutral atmosphere (factor ot 2 at auroral 
altitudes), the resultant chance in the modeled intensities is comparable to or less than the uncertainty in the 
neutral atmosphere. The smallest variations, for example, are for I 1838 (approximately 10 to 207c) while the 
largest variation is seen in the Ol 1356 A emission which is linear with 10] to within 20 7. We have also 
investigated the dependence of these intensities, and their ratios, to much larger changes in the composition (i.e. . 
[0]/|N-.]) such as might be encountered in large magnetic storms, or over seasonal or solar cycle extremes We 
find that the variation in ihe ( 1356/1 1838 ratio over the equivalent of a solar cycle is less than 507r. The 
summcMo-wmier changes arc approximately a factor of 2. The II 356/1 1 838 ratio is a very sensitive indicator ot 
the characteristic energy, showing a change of 1 3 over the energy range 200 eV to 10 keV The corresponding 
change in the LBH long-to-short wavelength ratio is much less (about a factor of 3). However, the latter is 
insensitive to changes in the neutral atmosphere ( <207 changes in LBH emission ratio for large changes in N : ). 
The three emissions therefore potentially provide a most valuable diagnostic of particle characteristic energy and 
energy flux. 


i . Introduction 

While in situ observations of energetic particles provide 
accurate information on the particle characteristics at the point of 
measurement, imaging from space ot the enure auroral oval holds 
the potential for providing details on total auroral energy influx, 
estimates of the characteristic energy of the auroral particles, and 
the capability to map and relate the footprint of this derived 
information hack along the magnetic field lines to various regions 
of the maenetosphere Auroral imaging in the vacuum ultraviolet 
permits observations of the regions ot interest under both day and 
nicht conditions. Work by Rees and Luckey [1974] on the ratios 
of visible emissions. UV emission intensity calculations by 
Strickland et al. [1983], and analysis of UV auroral spectra by 
Ishimoto et al. [ 1 988 1 al! indicate the potential value of using 
ratios of emission intensities to study auroral processes. A major 
focus of work in this area at the present time is to establish the 
quantitative footing on which such determinations can be placed. 

With the exception of HI Lya, the OI multipiets at 1304 A and 
1356 A and the N 2 Lyman Birge Hopfield (LBH) bands are the 
most prominent vacuum ultraviolet auroral emissions. The Ol 
1304 A emission has a high efficiency for multiple scattering. As 
a result, it has limited use for actual auroral imaging, although it 
does have potential value as an indicator of the O concentration. 
While the 1356 A emission does undergo multiple scattering, the 
efficiency is relatively small [Strickland and Anderson. 19831 and 
we ignore multiple scattering for I 1356 for this study. Similar 
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considerations allow us to also ignore multiple scattering for the 
Ni LBH emissions that are also considered in this study. The OI 
1356 A emission is absorbed increasingly by 0 2 with decreasing 
altitude. Thus its intensity vanes strongly (inversely) with 
increastne depth of penetration ot the incident auroral electrons 
and hence with increasing energy. The N : LBH transitions are 
electric dipole forbidden and the only prominent excitation mech- 
anism is electron impact. The LBH emission may therefore serve 
as a direct measure ot the total energy flux ot charged particles 
into the atmosphere. The longer wavelength LBH bands, which 
lie outside the region of substantial 0 2 absorption, are useful 
indicators of the total energy influx, while the long-to-short 
wavelength LBH intensity ratio provides information on the CL, 
and thus also some information on energy. These are the 
emissions (OI 1356, long and short wavelength LBH) on which 
we shall concentrate in this study. 

The purpose of this paper is to examine the sensitivity of these 
emissions to both likely uncertainties and anticipated changes in 
the neutral atmosphere. This is just one step in the process of mak- 
ing quantitative interpretations of auroral images, but an 
important one. We will consider other aspects (energy spectral 
characteristics and wavelength spectral extraction) elsewhere. In 
this paper we conduct a series of sensitivity studies using an 
auroral emission code that has been developed by our group 
[Richards and Torr. 1990]. The results are discussed below. 

2. Description of Auroral Code 

The behavior of auroral OI 1356 and N : LBH emissions has 
been studied with the use of an auroral computer model. The 
model is a two-stream auroral electron energy loss code that deter- 
mines the energy degradation of the primary spectrum as a func- 
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lion of energy and altitude and determines the production rates of 
prominent auroral emissions. A more complete description of the 
program, the selected cross sections, and comparison with other 
auroral models is given by Richar d and Torr (1990). The 
N : ( a 1 Try, ) cross section, as well as tl ' *auos of the individual LBH 
bands arc taken from Ajeilo and Sncmansky [1985]. Attenuation 
due to molecular oxygen absorption is explicitly computed with 
O absorption cross sections taken from Ogawa and Ogawa 
1 19751 and Hudson 1 1 971] The model currently assumes that 0 2 
absorption at wavelengths beyond 1750 A can be ignored. 

The model is optimized by incorporating variable energy bins 
\Torr ei al, 1974; Swartz. 19851 for the energy grid. To prevent 
numerical instabilities and violation of energy conservation, the 
altitude grid is variable to allow small grid steps (less than I km) 
at low altitudes. As a result, energy is typically conserved 
to within 10%. The code utilizes either the MSIS-86 neutral 
atmosphere \Hedin, 1987] or a user-supplied atmosphere. Either 
monoenergctic fluxes or a specified energy spectrum may be 
used. The incident energy spectrum may be modeled as a 
Gaussian or Maxwellian distribution, after Strickland et al. 
|19831 (hereafter SJ\V), or a user-supplied distribution may be 
used. All simulations reported below employed a I erg cm : s 1 
Gaussian incident flux distribution. The Gaussian scale parame- 
ter. labeled W in SJW , has been set equal to 0.25 E char which 
yields a full width at half maximum of 0.5 fln21 1 2 E char , where 
E char is the characteristic energy. 

3. Sensitivity Studies 

The emission studies reported here involved modeling auroral 
emissions at local midnight at 60 degrees north latitude. Table 1 
lists the MSIS parameters used in this study as well as detailing 
the range of solar activity investigated in the latter part of the 
studv 


TABLE ! MSIS Model Parameters 


Solar Activity 


Minimum 

Moderate 

Maximum 

1 ,,, - cm Flux Index 

75 

IK) 

200 

\\erage F, (1 - cm Index 

75 

J 10 

200 

\p Magnetic Activity Index 

4 

20 

1(H) 


Geographic latitude, 60 degrees: geographic longitude. 0 degrees: solar 
apparent lime. 0.0 hours: days 173. 356. 


Three emission ratios were studied. The first ratio was OI 1356/ 
LBH| (inL> , where LBH (oni . designates an N 2 LBH emission not 
strongly dominated by 0 : absorption. Specifically, the (2,8) band 
at 1838 A was chosen for this purpose. Second, the ratio O! 1 356/ 
LBH, hon was also modeled to investigate the relative influence of 
absorption by molecular oxygen. Here, LBH vhon is represented by 
the ( 1 , 1 ) band at 1464 A. By analogy with the previous definition 
LBH short is an LBH emission which is strongly absorbed by 0 2 . 
The final ratio studied was LBH| 00g ,LBK hon . The volume emis- 
sion rates integrated over altitude give the surface brightness or 
column intensity of the emission which we shall designate I 1 356, 
I 1464. and I 1838. These are the intensities that would be seen by 
a nadir viewing instrument from above the emission layer. 

3.1. Sensitivity to the Uncertainty of a Single Constituent 

The first question we chose to investigate was the dependence 
of the selected emissions and emission ratios on the uncertainties 
at any given time in our knowledge of the neutral atmosphere. We 
have assumed for this purpose that if we base our calculations on 
the MSIS-86 model atmosphere, the concentrations of O, 0 2 , or 
N 2 at auroral altitudes may on the average be uncertain by as much 


as a factor of 2. There will be occasions on which the uncertainty 
will exceed a factor of 2, but typically it will be less. This is 
similar to studies performed in 57W, but extends the investigation 
to study the dependence of auroral emissions to each of the major 
atmospheric constituents. In addition, in section 3.2 below, we 
further extend the study to include larger compositional variations 
due to seasonal and solar cyclic variations. 

The unperturbed, or reference, atmosphere is an MSIS model 
for high solar activity at summer solstice. Figure 1. Figure 2 



Fig. I . MSIS-86 reference atmosphere used as the standard case in this 
study (Day = 173. F m i - 700. Ap = 1 00). 



Energy (kev) 

Fig. 2. Nadir viewing column brightnesses tl 1356. 1 1464. and 1 1838) 
calculated using the reference atmosphere shown in Figure I. The 
diamonds here and in the remainder of the figures show the selected 
Gaussian characteristic energies for the incident electron energy 
distribution. 
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shows ihc computed I 1 356, I 1464, and 1 1 838 intensities for this 
model atmosphere as a function of energy over the range 200 eV 
to 10 keV as determined by our auroral code. In what follows we 
shall compare these results with those obtained when each of the 
atmospheric constituents (O, 0 2 . N 2 ) is. in turn, individually 
multiplied by 2 at ail altitudes, while the other tw r o are held 
constant. 

Before examining the results of these atmospheric changes, let 
us consider the possible impact of the changes. Doubling the con- 
centration of a constituent might at first be expected to double the 
effect of that specie on the column brightness of the monitored 
auroral emissions. In reality, however, there are a number of 
possible options in the interaction of the penetrating electrons and 
the atmospheric gases that render the situation more complex. An 
electron of a given initial energy will undergo a fixed number of 
collisions in a particular gas before themnaiizing. Thus in the very 
simple case of a single constituent atmosphere, changing the con- 
centration simply raises or lowers the altitude of the peak energy 
loss (and peak emission). This example (single-constituent 
atmosphere) is representative of those altitudes in which the con- 
centration of one atmospheric constituent dominates. From Figure 
! it can be seen that (for the conditions chosen) O tends to 
dominate above 400 km and N : tends to dominate below 300 km. 
In a mixture of gases, increasing the concentration of one specie 
relative to the others may also have the effect of raising the 
penetration depth. However, the gases will compete for incoming 
electrons in proportion to their mixing ratios and collision cross 
sections, and the ratio of the resulting emissions changes 
accordingly. 

Figure 3 shows the volume emission rate profiles for OI 1 356 A 
photons for the reference atmosphere case for selected energies. 
Only the very soft electrons (<200 eV) lose their energy in the 
altitude region where O is the major species. AM other energies lie 
in altitude regions where the various gases can compete for 
collisons with the precipitating electrons. Thus, for example, 
increasing the concentration of N 2 will result in a decrease in the 
production of OI 1356 A photons, because electrons that would 
have collided with O atoms now have an increased probability of 
colliding w'ilh N : molecules. The results of changing the O. 0 2 . 
and N : concentrations individually by a factor of 2 are shown in 
Figures 4 and 5. 



Production Rote (cm" 1 s* 1 ) 


Fig. 3. Volume emission rate profiles for OI 1356 A for ihc reference 
atmosphere shown in Figure 1 . 
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Figures 4</ and 5a show- the impact on the computed OI 1356 A 
and LBH surface brightness of the doubling of |N : ]. The 1356 A 
intensity at 2 keV drops to 65% of the reference model, and then 
rises back to 86% of the reference model by 10 keV. The I 1356 is 
reduced at all modeled energies due to the fact that there are now 
more collisions with N : . with a corresponding decrease in the 
production of O emissions. The decrease is not a full factor of 2 
because of the abundance of atomic oxygen at the higher altitudes. 
For the very lowest energies (highest altitudes), where (Nil is 
much less than |0], the 1 1356 should tend to an intensity level 
unchanged from the standard case, as the emission is simply 
raised in altitude. The modeled emissions do not include initial 
energies less than 200 eV which would lose their energy above 
400 km. but the lowest energy emissions do show this trend. At 
the higher energies, the emission is produced primarily at alti- 
tudes where N 2 is the major constituent and increased N 2 does not 
result in a significant change in the competition between O and 
NV 

For the LBH 1464 A emission, doubling the N 2 reduces the 
relative concentration of the dominant absorber. 0 2 . Thus for the 
higher energies which penetrate to greater depths, the emission 
from the increased N : overwhelms the 0 2 absorption. Absorption 
by 0 2 is not significant for the LBH 1 838 A emission; there is thus 
little dependence on the energy of the incident electrons. 

Doubling the 0 2 density (Figures 4 b and 5b) increases the 
absorption of I 1356 and I 1464 at the higher energies (lower 
altitudes) resulting in reduced column brightnesses. The LBH 
1838 A emission is relatively unaffected by 0 2 absorption and is 
influenced only by increased competition for collisions of the 
energetic particles with 0 2 molecules. However, since [N 2 ] 
remains the major specie relative to (0 2 ], I 1838 shows only small 
changes. 

The effect of doubling the [OI is show n for 1 1356 in Figure 4c 
and for the LBH emissions in Figure 5c. For I 1356 the effect is 
close to a factor of 2 increase in the emission at all energies, while 
for the LBH emissions there is almost no effect at ail energies. For 
the very low energies (not modeled) where O is the major con- 
stituent, the effect of doubling the O is simply to raise the altitude 
of the 1356 A emission. The trend to an unchanged emission can 
be seen at the lower energies. For the LBH emissions, the effect 
of doubling the O concentration is only seen at the very low 
energies (high altitudes) where the competition with N 2 is further 
increased. In the altitude regimes where N 2 is a larger component, 
the O does not play a significant role. 

From the results shown in Figures 4 and 5 it can be seen that 
uncertainties of a factor of 2 in any of the principal neutral atmos- 
pheric species translate into uncertainties ot less than 20% for 1 
1838. The LBH 1464 A emission shows variations up to 70% for 
factor of 2 uncertainties in [N 2 ] and less than 40% variation due to 
other constituents. OI 1356 is weakly sensitive to changes in 0 2 
and N 2 , but varies almost in direct proportion to changes in O. 


3.2. Sensitivity to Larger Compositional Changes 

The neutral atmosphere exhibits relatively large compositional 
changes in the course of the seasonal, solar cyclic and magnetic 
storm variations. In this section we report the results of our 
assessment of the dependence on the computed emissions chosen 
for this study to changes of this magnitude. In order to simulate 
changes we have varied the input parameters (F m7 , Ap, day of 
year) to the MSIS-86 model atmosphere, yielding neutral atmos- 
pheres at summer and winter solstice for conditions corresponding 
to low, moderate, and high solar activity (Table 1). The relative 
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Fig. 4. O? 1 356 dependence on (O] , [0 : ], and [N : ]. Each curve shows ( 1 356 for the doubled constituent case divided by T 1356 
from the standard (unchanged! case. 


compositional changes produced by these cases are illustrated in 
Figure 6. 

The variations in the emission ratio of OI 1356 A to LBH 1 838 
A due to such composition changes are illustrated in Figure 7. 
This particular ratio is very sensitive to incident energy. For any 
given atmospheric conditions, the ratio varies by a factor of 13 
over the energy range shown in the figure. What is interesting to 
note in Figure is that the variation due to the compositional 
changes produced by solar activity variations (low. moderate, and 
high Fjq 7 cm flux) are small (^30%), while the variations result- 
ing from compositional changes of the type produced by seasonal 
variations are much larger (about a factor of 2). In Figure 8 we 
show the energy dependence for these various composition cases 
for the individual I 1356 and I 1838 intensities. The LBH intensity 


is relatively insensitive to the changing atmospheric conditions, 
while the OI 1356 is found to be primarily responsible for the 
variations shown in Figure 7. 

The reason for this can be seen in Figure 9 which shows the 
ratio of the individual concentration changes relative to the 
standard case. Figure 9 a shows that the atomic oxygen concentra- 
tion (for altitudes below 300 km, which correspond to the initial 
energies modeled here) is significantly higher for the winter cases 
than the summer cases. 

Figure 10 shows the modeled volume emission rate altitude 
profiles for the solar minimum and solar maximum summer cases. 
As the atmosphere expands under the influence of increased solar 
activity, the production rate for a given energy peaks at higher 
altitudes. note is the fact that the behavior of the production 



Fig. 5. LBH dependence on [OI, [0 : ], and [N : ], Each curve shows the LBH intensity for the doubled constituent case divided by 
the LBH intensity from the standard (unchanged) case. 
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Fig. 6. Relative altitude density profile variations as a function of solar 
activity and seasonal variation. 



Energy (kev) 


Fig. 7. Sensitivity of the OI I356/LBH 1838 intensity ratio to solar 
activity and seasonal variation. Solid lines indicate summer conditions; 
dashed lines are for winter. Within each seasonal set of curves, the largest 
intensity ratios are obtained at solar minimum activity; the smallest ratios 
occur at solar maximum. 


rate profiles for electron energies below 2 keV is markedly dif- 
ferent from that above 2 keV. The altitude at which this energy 
loss peaks is approximately 140 km. This is the altitude below 
which 0 2 becomes a competitive constituent (see Figure 1). 

Figure 1 1 illustrates the effect of local O? absorption. As would 
be expected, for emissions lying outside the region of 0 2 absorp- 
tion (LBH 1838), the production rates are unchanged by local 0 2 
absorption. Foremissions within the Schumann-Runge absorption 
continuum, however, the shape of the emission rate profiles is 
changed significantly due to local 0 2 absorption. 

The final ratio modeled was LBH 1ong /LBH sh0 n (Figure 12). The 
ratio of LBH 1838 to LBH 1464 shows a dependence on the 
incident electron energy that varies only slightly with solar 
activity. As above, this variation can be explained by the relative 
densities of N 2 and 0 2 . The observed variability from solar 
minimum to solar maximum is due to changes in the 0 2 column 
density and hence in the 0 2 absorption. 

4. Discussion 

We have shown the intensity ratio OI 1356/LBH| 0ng to be a 
useful diagnostic for determining the characteristic energy of the 
auroral particles using LBHi ong to be LBH 1838. The I 1356/ 
1 1838 ratio is a very sensitive indicator of characteristic energy, 
changing by a factor of 13 or more over the range 200 eV to 10 
keV, but this ratio can vary by up to factors of 2 with changes in 
the neutral atmosphere. Almost all the change is due to variations 
in I 1356. In addition, the I 1838/1 1464 ratio shown in Figure 12 
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Fig. 8. Dependence oi LBH 1838 and Of 1 356 on solar activity and seasonal variation. The curves have the same interpretation as 
in Figure 7. 
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Fig. 9 Variation of altitude density profiles as a function of solar cyclic 
and seasonal variations. The profiles have been normalized to the standard 
case (summer, F, o 7 = 200). 


is found to be another potentially useful determinant of 
characteristic energy. The LBH long-to-short wavelength ratio 
varies by about a factor of 3 between 0.2 and 10 keV but, unlike 
the OI 1356 emission, is almost insensitive to changing atmos- 
pheric composition. The longer wavelength LBH bands are also 
useful indicators of the total energy influx, while the long-to-short 
wavelength LBH intensity ratio provides information on the O*. 

Other workers have also attempted to use optical emissions to 
characterize the auroral electron precipitatio ~or example, 
Ishimoto et al. [1988] studied a similar ratio (OI l_o6/LBH 1928) 
using a recent version of the Strickland model. Their modeled 
emission ratio shows a sensitivity of 9 between 1 and !0 keV, in 
good agreement with our results. In their description of the auroral 
code used in this study, Richards and Torr [1990] conduct a 
comparison between the two-stream model used here and the 
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Fig 10. Emission production rate altitude profiles tor Ol 1356 and LBH 


more sophisticated multi-stream model of SJW and find good 
agreement in the shape of the OI 1356 emission curve. (Dif- 
ferences in magnitude are due to the use of revised OI cross 
sections in our model.) In addition to these studies, Rees and 
Lummerzheim [1989] have suggested the N : 337 1 /N 2 " 4278 emis- 
sion ratio as a determinant of the incident auroral energy. Their 
results, however, disagree with calculations by SJW and with our 
model, which shows the 3371 emission to be independent of 
characteristic energy above 0.5 keV \ Richards and Torr , 1990]. 


We have investigated the sensitivity of OI 1356 A, LBH 1464 
A. and LBH 1838 A auroral emissions to changes in the neutral 
atmosphere. Our studies show that OI 1356 varies linearly with 
[O] to within 20 7c and shows much less variation with other 
atmospheric constituents. The LBH 1838 A intensity is relatively 
insensitive to typical uncertainties in the neutral atmosphere (fac- 
tor of 2 at auroral altitudes). 1 1464 shows larger variations 
because of its additional interaction with 0 2 . Our results are in 
good agreement with similar sensitivity studies performed by SJW 
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Fig. ! 1 . Same as Figure 10 (solar maximum), but illustrating the effects 
of local CK absorption The dashed curves are with local 0 : absorption. 


who used a Jacchia model atmosphere ( Jacchia, 19771 to model 
OI 1356 dependence on [01 The dependence of these intensities 
on much larger changes in the composition such as might be 
encountered over seasonal or solar cycle extremes has also been 
investigated. It is found that the OI 1356 A intensity is sensitive to 
compositional changes while tc N 2 LBH long wavelength emis- 
sion is relatively insensiti\ , such changes. 



Fig. 12. Sensitivity of LBH, imp /LBH, h<m ■ntjnsity ratio to solar activity 
and seasonal variation. The curves have the same interpretation as in 
Figure 7. 
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The capability has been developed to model thermospheric airgiow emissions on a semiglobal scale (L < 5). 
This model produces volume emission rates as a function of altitude, latitude, longitude, and local time tor any 
selected date, and solar and magnetic conditions. The model can thus be used to provide three-dimensional maps 
ot the selected emission lor comparison with data obtained from orbiting vehicles. As such it becomes an 
essential tool in the ptanmng and interpretation of airgiow observations. A unique leature ot the model is that it 
incorporates lull intcrhemisphenc coupling by solving all the appropriate coupled equations along the magnetic 
tlux tubes from the mesosphere in one hemisphere to ihe mesosphere in the other hemisphere. As a result the 
effects of conjugate photoelectrons (and heat Huxesi can be fully explored. In this paper we select two thermo- 
spheric emissions with which to demonstrate the capability. The first is the 7320-A emission from the metastable 
0~( : P1 The second is the permuted emission at 3371 A from the N : second poslive 0-0 band. These two 
emissions, for which the photochemistry is relatively well understood, are used to show the seasonal, diurnal and 
solar cyclic variations on a scale that covers mid- and low-latitudes, and the effects ot interhemispheric coupling 
1 conjugate phoioelectrons). 


Introduction 

Airgiow emissions arc important indicators of atmospheric 
composition and the mechanisms responsible for the production 
and loss of the particular excited state from which the airgiow is 
radiated. For example, the (>-0 band of the N : second positive 
system, which radiates at 3371 A. is excited in the airgiow by 
photoelectron impact, and is lost only by radiation. As a result, 
this emission is an excellent indicator ot the photoelectron excita- 
tion rate j Koppet al.. 1977], The 0~ l~P) nieiastable state, which 
radiates at 7320 A. is excited both by photoelcetrons and by 
extreme ultraviolet photons. However, because it is a long-lived 
state, it is lost by quenching by O. N : . and electrons in addition to 
radiation Thus at high altitudes, where radiation is the dominant 
loss mechanism, the 7320-A emission can be used to infer either 
the atomic oxygen concentration or the solar ultraviolet flux if the 
other is known \Meriwether et a/.. 1978; Rusch et a!.. 1976). 

In the past, a number of detailed studies of these emissions have 
been made using measurements of the surface brightness altitude 
profiles obtained from orbiting spacecraft (see. for example. 
Walker et ai. (19751 and Rusch et al. (19771). In this study we 
have globally modeled these two emissions for a variety of 
conditions. 

The model we have used here is one that we have steadily 
developed over the years. We solve the coupled time dependent 
energy, momentum, continuity, and photoeiectron transport 
equations from 80 km in one hemisphere, along a field line to 80 
km in the other hemisphere (Figure I) [Young et aL. 1980 a.b: 
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Richards and Torr . 1985a. 1 988 1 . The equations that are solved 
arc summarized below: 

!. Ion continuity equation for major ions is given by 


/W, 

df 


a-ZvVf-V < b, 


(I) 


where N, is the concentration of the nh major ion. Q, and L, are its 
production and loss frequency, respectively, and <i>, is the ion flux 
defined below. The electron density is assumed to be equal to the 
sum of the ion densities. 

2. Momentum equation is given by 

<J>, = N t U, (2) 

where 


V7-, + 

'•f 




(3) 


where U„ is the neutral wind velocity and where the subscripts i 
and j are applied First to 0 + and H* respectively* and then 
to*FT and He* respectively. 0*, hP , and He* arc coupled 
through collisions and Ihe polarization electric fields. Since the 
influence of O* and He* on H* becomes significant in dif- 
ferent altitude regimes and since the effect ot He * and O* is small 
due to their mass ratio, this decoupling of a system of three major 
ions into 2 pairs of major ions considerably simplifies the numer- 
ical calculation of the major ion densities. Here a,, and are 
thermal diffusion coefficients and D, is the ordinary ion diffusion 
coefficient of species /, as in the work by St. Maurice and Schunk 
II977|. 
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3. The thermal electron energy equation is 


and 


3 w L — = -NJTS U,-- N,W,.ST r -Vq, + lQ,-lL r 




ell 


4. The ion energy equation is 

3 N* El = - N,kT,V U, - ^ NJeU^T, -V-t/, + ~Z.Q, — 1L, 
2 d/ 2 

5. The ion heat flow equation is 


*" K 


f-x.vr,- ■ \ y vrA 

V w / 


(4) 


(5) 


(6) 


6. The electron heat flow equation is 

q f = -\ r V7V (7) 

The thermal conductivity coefficients X, and X,.. and the term 
5 = v,/v 

where v Jt ’ , v,/ , v,\ and v/ are the effective collision frequencies, 
are given by St. Maurice and Schunk |1977|. The ion-neutral 
collision frequencies v in are from Schunk and Nagy [1980|. 

7. The photoclcctron Liouville equation is given by 


s.t- 1 r - -7' : 4> t +7-,4»- + 

T,4>-+Ti<b* + 


ds B 


2<cos 0> <cos 0>’ 

<1 4 , 


ds B 


2<cos 0> <cos B> 


( 8 ) 

(9) 


where 


* + {E.s) 
$~{E,s) 
q{E.s) 


<cos <J>> 
B 


photoelcctron flux outward along s; 
photoeicctron flux inward along s\ 
photoelectron production rate in the range E to 
E + dE due to direct ionization processes: 
photoelcctron production in the range E to E + dE 
due to cascading from higher energy photoelcctrons 
undergoing inelastic collisions; 
average cosine of pitch angle; 
magnetic field strength; 


T 2 = s +/»/a,*]; 

k 



Fig. I . Illustration of (he interhemispheric nature of the code in which the 
coupled and time dependent equations are solved from 80 km in one hemi- 
sphere. along a field line, to 80 km in the conjugate hemisphere. 


n k Ath species number density; 

p t } photoelcctron backscattcr probability for elastic con- 
ditions with the Ath species; 

a/ photoclcctron total scattering cross section for elastic 
conditions with the Ath species; 

a,/ inelastic cross section for excitation of the Ath pan- 
icle species. 


The model includes an option to increase the 0"-0 collision 
frequency as recommended by Burnside et aL [ 1 987] . The use ot 
the Schunk and Nagy |I980| values here does not significantly 
influence the results of this paper. 

The above formulation corresponds to conditions where the dif- 
ferences between species temperatures and flow velocities are 
assumed to be small, i.c.. stress and nonlinear acceleration terms 
are neglected. In addition, density and temperature gradients 
normal to the geomagnetic field lines are neglected and we 
assume that the electron and ion temperature distributions are 
isotropic. 

The continuity equation is solved using a rather unique 
approach. We can rewrite ( 1 ) in terms ot a function F : 

JF = ^ + V ■ (N,U I )-Q, + L,N, (10) 

of 


Then, using a Newton iterative procedure to find the minimum of 
F, we solve for the density at the grid point, j. Figure 2 shows 
how the field line is divided into elements about the actual grid 
point j. The lower limit of the element (0 is placed midway 
between the grid points j and j- !, and the upper limit (w) is pla- 
ced midway between j and j + I The iower limit of one element is 
the upper limit of the preceding element. We then integrate (10) 
between u and 4, and obtain the densities by solving 




(ID 


where B is the amplitude of the magnetic field The values of rW/ 
dt % Q t and L arc obtained at the limits of integration, u and 1, by 
interpolation between the actual grid points. 

In the past, manv models have encountered numerical problems 
above about 3000 km. due to the large diffusion coefficient in this 
region which results in small density changes producing large 
chances in velocity. At lower altitudes, both ion-ion and ion- 
neutral collisions arc important, while at greater altitudes, colli- 
sions become less important and the plasma can be accurately 



Fig. 2. The numerical grid scheme used for (he solution ot the continuity 
equation. 
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described using a diffusive equilibrium approach. The solution of 
(II) distinguishes this model from previous methods which eval- 
uate the terms of the integral at only one point. As can be seen 
trom til) and Figure 2. the flux at the upper limit of one element 
becomes the flux at the lower limit of the next element. These 
fluxes must be identical and hence this method has been cal led the 
"flux preserving scheme.'' Furthermore, the flux at any grid point 
is closely tied to the flux at neighboring grid points, allowing 
stable solutions even for regions where large changes would give 
rise to unstable solutions with other numerical methods. 

The model starts at noon with "best guess" initial values. It is 
allowed to run for 12 hours (in local timet before results arc used 
in order to reduce dependence on initial conditions. It has been 
found that the ionospheric densities over a midnight to midnight 
diurnal cycle typically repeat with only small differences due. for 
example, to plasmaspheric refilling. The plasmaspheric H" and 
He " contents arc initially low. and the flux tubes are allowed to fill 
continuously. The transport equations for the three major ions. 
He". H". and O". arc carried out in two steps. As 
mentioned above, the coupled O" and H" equations arc solved 
first, followed by the He" and H" equations, where the latter 
use the O ‘ -H " results. Thus the equations for the three major 
ions arc essentially solved in a simultaneous manner Below 
approximately ISO km. NO" and 0>" become major ions, but 
arc obtained from photochemical equilibrium calculations. 

The numerical solution of these equations, boundary conditions 
and other details arc further discussed by Yotant et al. \ I980</./>|. 

The full interhcmisphcric coupling is difficult to handle, but 
once incorporated it imposes no artificial upper boundary con- 
ditions for both thermal and photoelectron fluxes. This is most 
important for thermal coupling and the proper treatment of con- 
jugate photoclcctrons [Richards and Torr. 1 985/? | which can be 
significant in the calculation of airglow' emissions. Typically, 
attenuation of the conjugate photoelectron flux by Coulomb colli- 
sions results in approximately 5 r r energy loss. The model 
includes the option to specify loss due to pitch angle scattering 
which we assumed to be zero for this paper. A tilted dipole 
approximation is used for the Earth s magnetic field [Richards 
and Torr. I986*/|. 

The concentrations of the major neutral species arc provided by 
the MSIS-86 | Hcdin. 1987| to the model which then computes the 
concentrations of minor and excited slate species and major ions 
| Torr. 1985: Richards et al , 1982*/. 1986/fl. In this paper a 
simple model giving daytime poleward wands and equatorward 
neutral winds was employed. This behavior is consistent with the 
results obtained using the method of Richards and Torr |I986/?] 
and Miller et al. [1986|. The model includes the option to use 
model winds of Hcdin et al. |1988|. Killeen er al. [ 1 987 1 and 
Killeen (private communication 1989). Note, if the Hcdin er al. 

1 1988 1 model is used, the O " -O collision frequency of Burnside et 
al. 1 1987| should be used in order to produce the observed h„,F2 at 
night. 

The chemistry of all significant emitting species is incorporated 
in detail, including the excitation of the metastabie states [M. R. 
Torr and Torr, 1982] and vibrational states {Richards et al.. 
1986*/: Richards and Torr, 1986c], and the odd nitrogen 
chemistry [Richards er at.. 1981, !982Z>: Richards, 1986]. The 
vibrational population distributions of N : are determined, an 
important factor in calculating the ionospheric O" and N t . 
concentrations. The model also includes calculation of ihe vibra- 
tional populations of N;>" , but this does not significantly affect the 
results reported here. The photochemistry' is that described by 


Torr 1 1985] as updated and shown here in Table I and illustrated 
in Figure 3. At equatorial latitudes, where electric fields play an 
important role, the electron concentrations arc obtained from the 
fully analytical ionospheric model of Anderson et at. [ 1 989 1 . 
Elsewhere, the electron densities arc computed self-consistemly 
by this model. The transition occurs between L - 1.5 and 1.8 
with interpolation between these L shells. The major elements of 
the code arc shown in Figures 4*/ and 4 h. The solar EUV flux is 
obtained in the following w r av. For solar minimum (FI 0.7 = 71) 
the model utilizes the F741 13 reference spectrum from Torr et at. 

{ 1979], with the fluxes below 250 A doubled as recommended by 
Richards and Torr |I984| and supported by Oga\va and Jud#e 
]1986|. For other levels of solar activity, each of the 37 
wavelength intervals is scaled linearly as a function of FI 0.7 
using the solar maximum measured flux at FI 0.7 = 206 given by 
Torr and Torr ( 1985], For further details, see Richards and Torr 
[19881. which also provides the cross sections used. 

The main outputs of the model include ion densities (O", 
O "< 4 M. 0"( : 0). O "( : F). H". He". N" . NO", N : ". 
Ni**). neutral densities (N( 4 5). N (~D). N {~P). NO. 0( [ D). 
0( l 5). NmA%"). NyM. electron and ion temperatures and 
flow velocities, the photoelectron flux, and a large number of 
emissions (sec Figure 4o>, 

During the Atmosphere Explorer C. D. and E program, 
numerous studies were conducted which compared the photo- 
chemistry of the code with in situ measurements of species con- 
centrations. The photochemistry' yields results consistent with the 
data base taken over the lifetime of the AE satellites. Generally, 
very good agreement with measurements has been obtained with 
regard to all parameters with the exception of high altitude 
electron temperatures in the plasmasphcrc and ionosphere |/Veir- 
herrx er at., 1989]. The model has been extensively tested against 
comprehensive satellite and incoherent scatter radar data bases 
j Yount? et *//.. 1980*/, b\ Chandler er al.. 1983: Richards and 
Torr. 1988: Scwherrx et al .. 1989: Horwirz er al., 1990: Richards 
ct al.. 19891. Apart from ihe input parameters (such as data and 
location i the only free parameter in ihe code is the pitch angle 
^altering of photoclcctrons in the plasmasphcrc. In its present 
form, (he model is ideally suited tor studies ot the airglow 
emissions. 

The code is run on the Marshall Space Flight Center CRAY 
XMP computer. Values of output parameters arc provided on a 
global grid of points, providing results in a latitude, longitude, 
altitude, and local time mesh for any selected date, or solar or 
magnetic conditions. For the cases discussed in this paper, wc 
have run the model for 144 flux tubes, which corresponds to 144 
northern and 144 southern hemisphere locations. Figure 5 shows 
the locations of the field lines along which the equations are 
solved. Those arc constrained to L ^ 5. 

We have chosen the 0"( : F) emission at 7320 A and the N : 
second positive 0-4) band emission at 3371 A for the initial global 
modeling. The calculations have been made for November 28. 
1983 for which the FI 0.7 cm flux w^as 89. and the Ap index was 
23. Thus the November 1983 calculations correspond to a period 
of relatively low solar activity. An earlier example of the results 
for the 7320 A case has been shown by Torr et at. \ 1990]. The 
results shown here represent a significant improvement over the 
Torr er at. |1990] case, in that wc have added approximately 50 
more flux tubes at low latitudes (144 versus 96). In addition we 
have used much smaller time steps through the twilight conditions 
(5 minutes versus 20 minutes). For comparison, wc have also run 
the calculations for the same day of year, but for conditions 
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TABLE 1. Summary of Photochemistry Used in the Interhemisphenc Model 


Reaction 

Number 

Reaction 

Rate Coefficient (cm 3 s' 1 ) 
or Rate (s _l ) 

Reference 

1. 

O " 4 + hv 

—4 x I0-' 2 (77300)° 7 

Torr (19851 

2, 

0 2 " -f- e~-+0 4 O 

1 .6 x ICT 7 (300/7,)° 77 for 7, 3= 1 200 K 
2 x I0- 7 (300/7,)° 7 for 7, < 1 200 K 

Torr and Torr (1981); 
Mehr and Biondi [ 1 969] 

3. 

O" 4 0 : -40 : " 4 O 

2.1x10-" {7„ + 2 7,/3 x 300}-“ 767 

Chen et al. [1978f 

4 

O" 4 N. + NO" + N 

1 .533 x I0"' 2 - 5.92 x I0- 1 ’ (7„/300) 
+ 8.60 x 10-' 4 (7/,/300) 2 

Sr. Maurice and Torr (19781; 
Albritton [1978]; 

ORiGsMAL PAGE is 

OF POOR QUALITY 

for 300 s 7, ff « 1700 K 
2.73 x I0"' 2 - 1. 155 x 10"' 2 (7„/300) 
+ I .483X I0 r,j (7 eff /300) 2 
for 1700 < T,„ < 6000 K 

Chen et al. (1978] 

5. 

NO" + e + N 4 O 

4.3 x ir 7 (7/300)-' 

Torr and Torr (1979] 

6. 

N 2 " + O+NO" 4 N 

1 .4 x ir'° (7,/300r 44 for 7, < 1 500 K 

McFarland et al. ( 1974]; Torr ( 1979] 

7. 

N 2 " + *--*N 4 N 

2.7 x 10- 7 

Abdou et al. ( 1984] 

8. 

N 2 " + O + O" 4- N 2 

0.07 * # (7/300)° 21 

McFarland et al. [1974] 

9. 

N 2 " + 0 2 +0 2 " 4 N z 

9.1 x I0-" exp (-0.002 T e „) 

Undinger et al. [ 1974] 

10. 

N" 4 0 2 + 0 2 " 4 N 

4 x 10-'° 

Huntress and Anicich [1976] 

II. 

N " + 0 : +N0" * O 

2 x I0 _l ° 

Huntress and Anicich ( 1976] 

12. 

N* + O.-t-NO* + O I'D) 

TS 

II 

O 

Langford et al. ( 1985] 

13. 

0 2 " + N -> NO " 4 o 

1 .2 x I0" l ° 

Fehsenfeld [1977] 

14. 

0 2 " + NO+NO" 4 0 2 

4.4 x 10-'° 

Lindinger et al. (1974] 

15. 

O *( 2 D) + N,-*.N 3 * + O 

8 x ICT 10 

Rowe et al. [ 1980]; 

Johnsen and Biondi ( 1 980] 

16. 

CT( 2 D) + 0-»-0 *< J S) + O 

5 x I0-' 2 

Abdou et al. [ 1984] 

17. 

0"( 2 D) 4 0 2 +0 : " 4 0 

7 x 10-'° 

Johnsen and Biondi ( 1980] 

18. 

0"( 2 D) 4 e '+0"( 4 S) 4 e ~ 

6.6 x I0-* (300/7,; 7 

Henry et al. (1969] 

19 

o*( 2 n -*-o*( ! d) + Aw 

0.173 s*' 

Seaton and Os ter brock (1957] 

20. 

0*( 2 P)-»-0*C‘S) + hv 

0.047 a - ' 

Seaton and Osterbrock ( 1957] 

21. 

O "( 2 P) 4 e‘ + 0"( 2 D) 4 t- 

i .5 x |(T 7 (300/7 r ) s 

Henry et al. ( 1969] 

22. 

0"( 2 P) 4 r-fO "( 2 D) 4 e' 

4.7 x I0~“ (300/7,) 7 

Henry' et al. (1969] 

23. 

0~( 2 P) 4 N 2 -4 products 

4.8 x IO -, ° 

Rusch et al. [ 1977] 

24. 

0"( 2 P) 4 O + products 

5.2 x I0-" 

Rusch et al. ( 1977] 

25. 

He 4 hv -4 He" 4 e~ 

4 Ox 10“* - 1 2x 10 7 (c) 

Torr and Torr [ 1985] 

26. 

He" 4 N,+N" 4 N 4 He 

i x icr 

Adams and Smith ( 1976] 

27. 

He" 4 N 2 + N 2 " ^ He 

6.5 X 10“’° 

Adams and Smith ( 1976] 

28. 

O" 4 H -4 H " 4 0 

2.2 x I0" n (T,) 5 

derived from Banks and Kockarts ( 1973] 

29. 

H" 4 0+0" 4 H 

2.5 x 10-" (7„) 7 

derived from Banks and Kockarts [1973] 

30. 

O" 4 N( 2 D)+N" 4 0 

5 x I0" M 

Torr et a!. (1979] 

31. 

NO" 4 e“ + N( : D) 4 0 

0 ■ ks where 0 = 0.76 

Kiev et al. [1977] 

32. 

N 2 " 4 *- + N( 2 D) 4 N 

0 ■ k 1 where 0 = 1.9 

Queffelec et al. [1985] 

33. 

Ni" 4 0 + N( 2 D) 4 NO" 

0 ■ kf, where 0=1.0 

Frederick and Rusch [ 1977] 

34. 

N" 4 0, + N( 2 D) 4 0 2 " 

0 • k w where 0=10 

assumed 

35. 

N( 2 D> 4 0+N( J S) 4 O 

~7x I0-' 7 

Richards et al. [1981] 

36. 

N( 2 D) 4 0> + N0 4 O 

6 x i cr 1 2 

Lin and Kaufman (1971] 

37. 

N( 2 D) 4 4 e~ 

5 x lO -10 (7/300) 5 

Frederick and Rusch ( 1977] 

38. 

N ( 2 D) 4 0 2 " + NO" 4 0 

1 x ICT 1 

Dalgarno [1970] 

39. 

N( 4 S) 4 0,-4 NO 4 O 

4.4 x lO' 12 exp (-3220/7) 

Becker et al. [1969] 

40. 

N 4 NO+N : 4 0 

3.4 x I0-" 

Lee et al. [ 1978] 

41. 

O, 4 hv - 4 CH f D) 4 0 

0 = 1;J X (0>) SR = (15- 2.8) x 10^ 

Torr et at. [1980] 

42. 

O," 4 e'-rOi'D) 4 O 

0 k 2 where 0 = 1.2 

Abreu et al. [ 1 986] 

43. 

0 CD) 4 N : + 0 <V>) 4 Nj 

2. Ox 10 _l 1 exp ( 107.8 T„) 

Streit et al. [1976] 

44. 

O('D) 4 0 2 + 0 <V>) 4 0 2 

2.9 x |(T n exp (67. 5/7, f ) 

Streitetal. [1976] 

45. 

0('D)+0( 3 P) 4 hv 

0.00934 

Fischer and Saha ( 1983] 

46. 

N i 2 D) 4 O 2 + O(‘0) 4 NO 

—5 x I0 -12 

Rusch etai. (1978]; D. G. Torretal. [1981] 

47. 

0 2 " 4 ?>0('S) 4 0 

0 k 2 where 0 = 0.08 

Bates and Zipf (198!]; Abreu et al. [1986] 

48. 

0 2 " 4 N+O ( l S) 4 NO" 

~2 x lO -1 1 

Frederick et at. [1976] 

49. 

O('S) 4 0( 3 P)+0 4 O 

2 x i<r‘ 4 

Slanger and Black [1981] 

50. 

O i'D) 4 0-40 4 O 

8 x ICT 12 

Abreu et al. [1986] 
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TABLE 1. (continued) 


Reaction 

Number 

Reaction 


Rate Coefficient (cm 3 s* 1 ) 
or Rate (s _l ) 

Reference 

51 

0< l S) -+CM 'D +■ hv 

1.07 


Kernahan and Pan% (1975] 


-+Oi*P) + hv 

0.0444 



52. 

N CD) + NO -> N : -t- O 

7x J0-" 

i 

Lin and Kaufman [1971] 

53. 

Of 5) + O z + OCP) + 0> 

4.9 x 10 

exp ( — 1 730/«r> 

Zipf [1979] 

54. 

+■ 0^- products 

2 x I0‘" 


Piper [1982) 

55. 

N : (A 3 I,r> + O+Ot'S) + N, 

[3 • k *4 where (3 = 0.37 

Piper [1982] 

56. 

N^* + 0+0^ + N> 

~2 x 10 

- 1 LM» 

Abdou et al. [19841; Torr [1985] 

57. 

N, + hv + + N + e 

(1.78 - 

5. 14) x IQ-'* 

Torr and Torr [ 1 985]' 



(3.06 - 

8.82) x )0' 7 


58. 

0 + hv-^O* (*S) + e 

(0.98 - 

2 . 8 1 ) x 1 0~ 7 

Torr and Torr [ 1 985 )‘ 


■> 0~( 2 D) + e 

(0.79 - 

2.34) x 10“ 7 



-+0 ~( 2 P) + e 

(0.45 - 

1 .38) x I0 -7 



-►OVP) + e 

(1.04 - 

3.43) x |(T* 



•+0~t 3 r) + e 

(0.46 - 

i ,42) x 10** 



"Since the results of Chen et al. {1978] stop at 700°K, we normalize the convened drift tube data parameterized by St. Maurice and Torr (19781 at this 
(emperature. 

'The model computes this rate coefficient as a function of Ni * vibrational temperature lT t ). The value reduces to that given by equation (6) when T , = T n . 
Inclusion of this process is an option available [see Abdou et al., 1984], 

'The ranges given for the ionization frequencies indicate the variation over a solar cycle. 



Fig. 3. Schematic of the thermospheric and ionospheric chemistry that is incorporated in the model. 
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corresponding to solar maximum if'10.7 = 194, Ap - 23). This 
study deals only with the airglow. and we do not show results tor 
latitudes higher than L = 5. 

Details of the results are given in the following section. 

Modeling ok the 7320-A Airglow Emission 

The 7320-A emission arises from the O state, which is 
produced in the thermosphere by tw'o mechanisms: photoioniza- 
tion and photoelectron ionization: 

O A hv -*■ O "( : P) + c 
O A e* ->0 ~( 2 P) A 2c 

where the photoelectrons may be those produced locally, or those 
produced in the conjugate hemisphere and transported along the 
magnetic field lines. The 0^( : P) is lost by several mechanisms: 

- idiation 

O (*P)-4 O ( 4 5) A hv-ni o 
collisional deactivation 

OVj P) A a O 

o~( : P) + o -^cr - o 

O ~{ : P) A <* O ~( : D) 4 e 

O'C'P) + v O ~( 4 5) A c 


This photochemistry has been reviewed by Ton ami Torr 1 1982|. 
It should be noted that since publication of the acronomicaliv 
derived values of the rate coefficients by Ritsch er al. [1977], the 
solar EUV flux below ^250 A was doubled. We estimate that 
when quenching dominates, the reported intensities may be about 
30# too high. 

Plate l shows the 7320-A volume emission rate at the peak of 
the lavcr as a function of latitude, local time (longitude) and alti- 
tude. The upper plot shows the peak volume emission rate as a 
function of latitude and local time. Because the model has long- 
itudinal variability resulting from both the MSIS | Heclin 1987] 
model atmosphere, and the interhemispheric coupling, this par- 
ticular plot is shown for (X) UT. This UT is equivalent to placing 
midnight at 0° longitude, the Greenwich meridian, and noon at 
180° longitude, i.c. over the Pac i tie Ocean. Where the volume 
emission rate becomes so small as to be effectively zero, the 
values arc not plotted. These regions can be seen near midnight 
for equatorial latitudes, and represent no production at these 
times. The graphics tend to smear the northern and southern lati- 
tude boundaries by a few degrees to the north and south, respec- 
tively. The solutions are only valid, however, within the L < 5 
region illustrated in Figure 5. 

The ha :c features shown in Plate 1 have been discussed by 
Ton ct al. |I990| but we will review- them here briefly as this 
plate w ill be used for the comparison with other cases. Summer is 



Fig. 4<oj: Schematic illustration of the major inpui and output elements of the field line interhemispheric plasma (FLIP) code, ib) 

Illustration of the flow of the solutions. 
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Fig. 4. (continued) 
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in i he southern hemisphere .o immcdiutcK seen bv the longer 
duration ot (he middav peak *...|ume emission rates. Figure b 
'Iiihu the prtKiueium and loss rates at northern and southern mid- 
latitudes near noon In both ^..ws the major production mech- 
anism is pnotoiom/aiion with nhotocleciron impact contributing 



0 4 8 ’2 16 20 

Local Tim* 


I ie s 1 -valioiis >M ji kl potnis mch dturnai alniude solutions were 
•hi.nned 


about 25*1 The maior high aintude loss mechanism is radiation, 
while quenching b\ N dominates at lower altitudes It is interest- 
ing to note in Plate I that the middav peak volume emission rates 
m the winter i northern i hemisphere are somewhat larger than in 
the summer hemisphere. Hie reason tor this can he seen in Figure 

0. W hile the \ peak quenching rale is almost identical in both 
hemispheres, the V quench ins: tails oil with altitude with a 
smaller scale height in the winter hemisphere, so that radiation 
dominates to a lower altitude with the resulting higher emission 
rate . 

At and bevond the terminators it'or solar zenith angles larger 
than 00M sunlight illuminates increasingly higher altitudes, and 
the peak altitude rises while the volume emission rate tails 
sharplv In Figure 7 we show the production and loss rates for a 
solar zenith ancle of 105° in the evening southern hemisphere. 
The peak L V photoioni/ation has risen to above 600 km. A 
second peak is formed near 300 km. The latter is a result of photo- 
electrons which are produced in the ioni/ation process that caused 
the upper peak. The photoelectrons are transported downward and 
lose their enemy near 300 km in the impact ioni7ution of atomic 
oxygen. Plate 2 shows the same information as was given in Plate 

1. but versus solar zenith angle instead of local time. This plate 
shows the extent to which seasonal hemispheric asymmetries arc 
reduced when the data are plotted m this format. 


7320 A Emission Rate (F10. 7 = 09.0) 



- 1.5 1.2 

-1 .5 1.5 

BELOW —1.3 


Pl.iie 1 . The 7320- A \ ulume emission rate i photons cm" s“’ ) at the peak ol the layer as a f unction oi altitude, latitude, and local time 
i longitude) lor (Hi FT The upper plot shows ihe peak volume emission rate only as a function of latitude and local time. The volume 
emission rate is on a log scale The results are only valid wiihtn the Z. - 5 boundaries shown in Figure 3 
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Fig. 6. (continued) 


An interesting feature of Plate 1 is the fact that there is signifi- 
cant production (~0.2 cm" 3 s' 1 ) at high altitudes at midnight in 
the winter hemisphere. This is due entirely to photoelectrons 
transported from the summer hemisphere. This conjugate electron 
production rate together with the loss processes are shown in 


Figure 8. Without the intcrhcmispheric capability this production 
source would not be modeled. 

Because of the longitudinal asymmetries of the Earth's mag- 
netic field, the conjugate photoelcciron production rate varies for 
different maps generated for different universal times. Plate 3 
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Plate 3. The 7320- A volume emission rate (photons cnf ' s -1 ) at the peak of the layer as a function of altitude, latitude, and local lime 
lor 06 UT The volume emission rate is on a log scale. The missing areas near the equator at midnight are where the emission rate has 
L’lTociivelv stone lo zero. 
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Fig. 7. Production and loss rates lor 54°S lor a solar zenith angle of 105° (LT = (X)}. 


shows the global plot for the same conditions as Plate I. but for are present all through the night near 50°N. Thus the conjugate 
UT = 0600. This is equivalent to placing midnight at 90°W. over source is seen to be significantly more pronounced than for 00 

the eastern United States. In this longitudinal sector, higher UT. Plate 4 shows the integrated volume emission rate, or surface 

southern latitudes (solar illuminated) map to lower northern lati- brightness lin raylcighs). with and without the conjugate photo- 

tudes. Above 50°S the Sun docs not set and conjugate electrons electron source. The scale has been adjusted tv nhancc the winter 
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Plate 4</ With conjugate photoelectron production. 



Plate 4 h Without conjugate photoclcciron production 


PLite 4 Vertical 7320- intensity irayleighsi lor solar minimum at IT - 06. The solution is only valid within the L ? boundaries 
shown in It cure > 
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Plate 5. The 7320- A volume emission rate (photons cm" ' s _l ) at the peak ot the layer as a function ol altitude, latitude, and local time 
lor (X) IT The results are lor (he same dav as Plate I ( November 28) but tor solar maximum conditions. The volume emission rate is 
■>himn on a log scale 


nighttime conjugate production. This mtcrhemisphenc source 
results in 2 to 3 R of 7320-A airglow. 

Plate 5 shows the same day of year as Plate 1 . but for conditions 
representing solar maximum. The principal differences between 
these two arc in the overall increase in emission rate and the 
overall increase in altitude of the emission peak. Both of these 
effects arc to be expected on the basis of the solar flux and neutral 
atmosphere changes. 

Plate 6 shows the results of integrating the volume emission 
rates of the solar minimum and maximum cases (Plates i and 5) 
over altitude to obtain the vertical surface brightness. These 
intensities tin ravleighs) arc what an orbiting instrument would 
sec looking vertically downward on the atmosphere. 

Modeling of the N> 337 1 -A Airglow 

The production mechanism for the excitation of the N : second 
positive system in the daygiow is photoelectron impact: 

N : (,V'1,T) - e* N : (Cy„! - c 

The 337 1 -A emission arises from the 0-0 hand of the second posi- 
tive svstom: 

N : (CV,) ) * hv 


Since the transition is permitted, radiation is the only loss process. 
The excitation cross sections for the 0-0 band are those ot hnami 
and Borst [1974| which are 25 *7c of the total cross section. 

Plate 7 shows the results of the semiglobal solution of the 3371- 
A volume emission rates. The production rates tor noon at mid- 
iatitudes arc shown in Figure 9 The results shown in Plate 7 
follow the behavior anticipated from a simple photoclcctron 
source and radiative loss. The behavior changes after sunset in 
that the peak production rises in altitude and tails in magnitude. 
Again, the high latitude conjugate photoclcctron source can be 
scon during the night in the winter hemisphere. While we do not 
show the 06 UT case here, the conjugate source again becomes 
more pronounced for the different magnetic field orientation, as it 
did for the 7320-A case, because of the more favorable alignment 
of the magnetic field lines tor this purpose. 

Plate 8 shows the global model of the peak 3371 -A volume 
emission rate for solar maximum. What is interesting to note in 
this case is that, apart from an overall increase in the altitude ot 
the emission peak, the peak emission rate distribution is very 
similar at solar maximum and solar minimum. However, the dif- 
ference in the vertical intensity is significant. This is seen in Plate 
9 which shows the integrated column surface brightness for solar 
minimum and maximum. The difference is largely due to the dif- 
fering scale heights, as can be seen in Figure 9. 
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Plate bb Solar maximum case shown in Plate 5. 

Plate 6. Vertical column integrated surface brightness tin ravleighsl. The results arc only valid w ithin the L ^ > boundaries shown in 

Fieure 5 
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PI ale 7 Volume emission rate (photons cm"’ s“ l ) at the peak ot the 3371 - A layer versus aimude. latitude, and local time ( longitude > 
tor near solar minimum conditions The volume emission rate is plotted on a log scale The results are onlv valid vuihin the L ^3 
boundaries shown in Fieurc 5 


Summary 

In this paper we have presented the first mterhcmisphcricaliy 
coupled models of the 0“( : P) and N : t'C) airgiow emission at 
7320 A and 3371 A. respectively, covering mid and low latitudes. 
These are two emissions for which the mechanisms arc relatively 
well understood, allowing us to examine here the momnoiogy. 
The senuglobai graphical maps clearly illustrate the diurnal, 
seasonal and solar cyclic variability of these emissions In both 
cases, conjugate photoelectrons result in excitation at mid to nigh 
nighttime latitudes in the winter hemisphere. The in 1 mi- 

spheric coupling capability of the model allows us to he 


extent and the longitudinal variability of the conjugate photoelec- 
tron effects. 

This '‘-odd provides a capability which will greatly enhance the 
ability tu interpret airgiow observations made from spacecraft. 
Any line of sight geometry can now he projected through the 
three-dimensional model solution, and by integrating along the 
projected line of sight, surface brightness (vertical or slant path) 
can be obtained for comparison with the observations. Further- 
more. because of the complex variability exhibited by these 
emissions over the L ^ 5 range, this model also provides a 
valuable tool for planning of such observations. 


Plate ^ Vertical intensity ot ihe 3371 -A airgiow (in rayteighs) for solar minimum and solar maximum conditions. The solar 
maximum plot is shown on two scales; one to show the intensity vanatton and one to allow comparison with the solar minimum plot 
shown in Plate Va The results arc only valid within the .* 5 boundaries shown in Figure 5. 
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Plate 8. Volume emission rate (photons cm’ 1 s _l ) at the peak ol the 337 1 -A layer lor solar maximum conditions. The volume 
emission rate is plotted on a log scale 
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Plate 9b. Solar maximum, scaled to brightest intensities. 



Plate 9r. Solar maximum, same scale as solar minimum. 
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A MIDLATITUDE INTERHEMISPHERIC MODEL OF THE 0*( ’P) AIRGLOW EMISSION AT 7320 A 

Marsha R. Ton- 

Space Science Laboratory, Marshall Space Fight Center. Huntsville. Alabama 
D. G. Ton* and P. G. Richards 

Center for Space Plasma and Aeronomy, The University of Alabama in Huntsville 


Abstract. The results are reported of the first interhemispheric 
model of the 7320 A airglow covering mid- and low latitudes, at ail 
longitudes. A comprehensive model of the ionosphere and thermo- 
sphere is used to compute volume emission rates as a function of 
altitude, latitude, longitude, and local time. Selected results are 
shown here to illustrate the computational capability. In particular 
we discuss the diurnal and seasonal variability and interhemispheric 
coupling of conjugate photoelectrons. The model is particularly well 
suited for airglow studies, and provides a valuable tool for the 
comparison, interpretation, and planning of spectroscopic observa- 
tions made from orbiting platforms for all significant thermospheric 
emissions. 

Introduction 

The O * (~P) airglow emission at 7320 A provides a means of deter- 
mining the O * f 2 P) concentration. Measurements of this emission can 
be used to determine the atomic oxygen concentration or the solar 
UV ionization frequency [Meriwether et al., 1978; Rusch et a!., 
1976]. The photochemistry of this species has been established 
largely on the basis of the comparison of steady-state altitude profile 
calculations with surface brightness measurements made from satel- 
lites [Walker et al. , 1975; Rusch et al. . 1977]. Observations of limb 
brightness are inverted to yield altitude profiles of volume emission 
rates. The principal sources and sinks of O ^ ( 2 P) have been discussed 
in detail by Torr and Torr [1982]. 

Over the years, we have developed a comprehensive model of the 
ionosphere, thermosphere, and plasmasphere. This model solves the 
coupled time-dependent energy, momentum, continuity, and photo- 
electron transport equations from 80 km in one hemisphere, along 
the field line, to 80 km in the conjugate hemisphere [Young et al., 
1980]. The full interhemispheric capability allows for the proper 
treatment of thermal coupling [Richards and Torr, 1986] and con- 
jugate photoelectrons [Richards and Torr, 1985], which is important 
in the calculation of airglow emissions. The concentrations of the 
major neutral species are provided as input from MSIS-86 [Hedin, 
1987] to the model which then provides the minor and excited state 
species and ions [Torr, 1985], The chemistry of all the emitting 
species is comprehensively included [Torr, 1985]. Ab initio calcu- 
lations of the excitation and loss rates are performed for the metas- 
table species [Torr and Torr, 1982] and the vibrational states of 
molecules and ions [Richards et al., 1986], The concentrations of 
odd nitrogen species are also computed together with those of other 
minor constituents [Torr et al . , 1 980] . The photochemistry currently 
used in the model is that defined by Torr [1985] (Tables A-8 through 
A- 1 2 and A- 19 through A-21) with the updates, corrections, and 
additions summarized in Table i . 

The three-dimensional capability is achieved by running the code 
for approximately 100 magnetic flux tubes. Because the model is 
interhemispheric, solutions are obtained simultaneously for the con- 
jugate hemisphere. This yields a total of approximately 200 grid 
points. The flux tubes are selected along various L-shells for L ^ 5 
and the grid on w-hich the code is run is shown in F/igure 1 in geo- 
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graphic space. Because of the large computational requirements, the 
MSFC CRAY II computer is used. Results are output on a four- 
dimensional global grid, comprising latitude, longitude, altitude, 
and local time. Input parameters are date and solar and magnetic 
indices. 



0 60 120 180 240 300 

LONGITUDE fE) 

Fig. I . Location of the conjugate points along L = I 08 to L = 5 
which mark the held lines along w r hich the computations were made. 


In this paperwe have selected the 7320 A emission to illustrate the 
global airglow modeling capability. The calculations have been 
made for conditions corresponding to November 28, 1983, for which 
the F10.7 cm flux was 89 and the Ap index was 23. Thus these calcu- 
lations correspond to a period of relatively low solar activity. 

Results 

Figure 2 shows the results of the global solution ot the 7320 A 
emission. The volume emission rate at the peak of the layer is plotted 
as a function of altitude, latitude, and local time (longitude). This 
particular plot is shown for 00 UT, which is equivalent to placing 
midnight at 0° longitude, and noon at 180° longitude. 

Several interesting features appear in the results shown in Figure 
2. Summer is in the southern hemisphere and the more extensive 
solar illumination in this hemisphere is immediately evident in the 
longer duration of the daytime peak values (volume emission rates of 
the order of 10 photons-cm ^-s' 1 ). The model results are valid for the 
latitudinal regime indicated by the envelope of the points shown in 
Figure 1 . The interpolation routine tends to smear the high latitude 
boundary in Figure 2 by about 3°. Figure 3 shows the production and 
loss rate profiles for noon at southern midlatitudes. The peak produc- 
tion of 0 + ( 2 P) occurs near 170 km with photoionization the major 
source. The dominant loss mechanism above 280 km is radiation, 
and below 280 km it is collisional deactivation by N 2 . The combina- 
tion of these processes results in the emission peak being formed at 
approximately 260 km. 
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TABLE I Updates to Photochemica! Model of Torr ( I9S5| 


Reaction 

Rate Coefficient (cm 3 s' 1 } or s' 1 ) 

Reference 

1) O' - o : ->o; - u 

. .0- .' T - - IT '| — 
( 3. 300 \ 

Chen et al. [1978] 

2l N. ' - O -»0; - N; 

9 1 x 10 " exp <-0 002T,.„) 

Lmdinger et al. [ 1974] 

3) O' - WDi *+ N ‘ - O 

5 x 10" 

Torr et al. [ 1979] 

4iO* - H -H * + O 

2.2 x 10 " T, 5 

Banks and Kockarts [1973] 

5) H * - O-CT - H 

2.5 x 10 " T„'5 

Banks and Kockarts [1973] 

Pa , 

6) NT ^ e- -► N( : D) + n 

a = 2.7 x 10' 7 
(3 = 19 

Abdou et al. [1984] 
Queffelec et al. [1985] 

7) O('D) + O; -Ol'P) + 0 2 

2.9 x 10 “e 675 T " 

Streit et al. [1 976] 

X) Ol 1 D> -OCP) - hv 

\, o = 0.00934 

Fischer and Saha [1983] 

9) N : (A^) +■ O -+ products 
-Ol'S) + N : 

2 x I0 M 
(3 = 0.37 

Piper et al. [ 1981 ] 
Piper [1982] 

10) Ni : D) - NO -*> N : + O 

7 x 10" 

Lin and Kaufman [1971] 

ll)0('D)* 040 + 0 

8 x I O' ' : 

Abreu et al. [ 1986] 

12) 0( S) -OLD) - hv 
OCP) 4 . h v 

1 .07 
0.0444 

Kemahan and Pang [1975] 


Fig. 2. Three-dimensional plot of the 7320 A volume emission rate at the peak ot the layer versus latitude and local time, tor 00 
UT The calculations were done for conditions appropriate for November 28. 1983. The upper panel shows the same results, 
hut without the altitude coordinate, i.e.. peak volume emission rate versus latitude and local time. The color code is on a log 
scale, and the units are photonscm ’ s ' 
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At twilight the situation changes rapidly. The photoionization 
peak rises rapidly in altitude and the production rate falls sharply. 
Figure 4 shows the production rates for 48°S at 22:23 LT. At this 
time the solar zenith angle is 107°, so that the high altitudes are 
illuminated. There is a second peak at about 300 km due to trans- 
ported local photoelectrons which are produced during the ionization 
process which created the upper peak. Radiative loss dominates most 
of the profile. At low latitudes at midnight, the peak production rate 
tends to zero at very high altitudes. When the production rate effec- 
tively drops to zero these values are omitted in Figure 2, so as not to 
obscure the main results. 
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In the winter (northern) hemisphere the twilights draw closer 
together but the photochemistry of the daytime and twilight 7320 A 
emission is essentially the same as that described above. What is 
interesting to note here is that there is some proo iction of 0"( : P) at 
high latitudes even at midnight. The source is impact ionization by 
photoelectrons from the sunlit conjugate hemisphere. The conjugate 
photoelectron source produces a column integrated surface bright- 
ness of 2 to 3 R of 7320 A emission. 

Figure 5 shows the global map of the 7320 A intensity for the same 
case as Figure 2. This is what would be seen by an instrument look- 
ing straight down at the Earth for a fixed UT. The midday intensities 
in the winter hemisphere are slightly larger than those in the summer 
hemisphere. This is because the winter hemisphere is somewhat 
colder, and radiative loss dominates to a lower altitude. 


Summary 

We have reported the first results of a midlatitude interhemispher- 
! v coupled model of the 7320 A airglow. These results show the 
i. .rnal and seasonal variation for conditions representative of near 
solar minimum. The magnitude and extent of the conjugate photo- 
electron source is quantified, together with the sources and sinks at 
selected times. In another paper [Torr et a!., 1989] we show further 
details of the longitudinal and solar cycle variations. This model 
allows three-dimensional maps of any airglow emission to be gener- 
ated. Any line of sight geometry can be projected through these maps 
and the volume emission rate integrated along the viewing direction 
u table comparisons with observation. A three-dimensional 
perspective of airglow emission rates or brightnesses is of relevance 
to the interpretation of measurements taken from an orbiting vehicle. 
A satellite flying in a circular orbit through the topology repre- 
sented in Figure 2 would observe large variations in emission 
through the structured regions. Semi-global maps of the type pro- 
duced here provide rapid insight into the sources of the variability. 
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NEW HIGH PERFORMANCE FILTERS FOR VACUUM 
ULTRAVIOLET APPLICATIONS 
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AND 
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OBJECTIVES 


TO DEVELOP NARROWBAND HIGH "TRANSMITTANCE" FILTERS FOR 
THE VACUUM ULTRAVIOLET (VUV) WAVELENGTH REGIME. 

GOALS: 

TRANSMITTANCE > 30% FOR 130 < X < 200 nm 
BANDPASS ~ 3 nm FWHM 
OUT-OF-BAND BLOCKING < 10 3 % 


APPLICATIONS: 

SPECTRAL ISOLATION OF THE TERRESTRIAL 
EMISSION FEATURES: 

HLya 121.6 nm 

01 130.4 nm 

Ol 135.6 nm 


AND OTHER FEATURES 


CONCEPT FOR THE DEVELOPMENT OF HIGH 
PERFORMANCE FILTERS 


PROBLEMS: 


ALTERNATIVES: 


PRINCIPLE: 


LACK OF LOW ABSORPTION FILM MATERIALS IN THE 
VUV FOR ALL-DIELECTRIC TRANSMISSION FILTERS 

POOR TRANSMISSION, BANDWIDTH AND BLOCKING 

REFLECTIVE COATINGS ARE NOT AS SEVERELY 
AFFECTED BY ABSORPTION EFFECTS 

BUT TYPICAL REFLECTIVE COATINGS HAVE A 
RELATIVELY BROAD PASSBAND 

A MEASURE OF SUCCESS HAS BEEN ACHIEVED 
PREVIOUSLY BY COMBINING SEVERAL REaECTORS 
IN TANDEM 

THE RATIO OF IN-BAND TO OUT-OF-BAND 
REFLECTIVITY (R) INCREASES AS 

(R in/Rout) N WHERE N = NUMBER OF REFLECTORS 


IMPROVED BANDWIDTH AND OUT-OF-BAND BLOCKING 


Our Objective and Commercially 
Available VUV Narrowband Filters 



X ( nm ) 


Peak Transmission 
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RGURE 1 : THE MAXIMUM REFLECTANCE OF THE n STACK CALCULATED FOR ZERO ANGLE OF 
INCIDENCE AT 145 NM. DIAMONDS REPRESENT 99-LAYER, TRIANGLES 55-LAYER, AND SQUARES 35- 
LAYER STACKS. H = LANTHANUM FLUORIDE, AND L = MAGNESIUM FLUORIDE. THE KOPPELMANN LIMIT 
IS 90.8%. 
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FIGURE 2: THE MAXIMUM REFLECTANCE OF THE n STACK CALCULATED FOR 45 DEGREES ANGLE OF 

INCIDENCE AT 135.6 NM. DIAMONDS REPRESENT 99-LAYER, TRIANGLES 55-LAYER, AND SQUARES 
35-LAYER STACKS. H = LANTHANUM FLUORIDE, AND L = MAGNESIUM FLUORIDE. THE KOPPELMANN 
LIMIT IS 89.5% 
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FIGURE 5: THE MEASURED (SQUARES) AND CALCULATED (DIAMONDS) REFLECTANCE OF THE 29- 
LAYER n STACK FOR 45 DEGREES ANGLE OF INCIDENCE CENTERED AT 135.6 NM. THE OPTICAL 
THICKNESS RATIO H/L - 1/3 WITH H = BARIUM FLUORIDE, AND L = MAGNESIUM FLUORIDE. 
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RGURE 6: THE MEASURED (SQUARES) AND CALCULATED (DIAMONDS) REFLECTANCE OF THE 35- 
LAYER SECOND ORDER QW STACK FOR 45 DEGREES ANGLE OF INCIDENCE CENTERED AT 1 35 6 NM H = 
BARIUM FLUORIDE AND L = MAGNESIUM FLUORIDE. 
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RGURE 7: THE TRANSMITTANCE OF THE COMBINATION OF FOUR 29-LAYER FILTERS SHOWN 
IN FIGURE 5. THE BANDWIDTH IS 4.3 NM, AND A PEAK TRANSMITTANCE AT 135.6 NM IS 53.7%. 
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Multiple Reflectors as 


Narrowband and Broadband VUV Filters 


Muamer Zukic and Douglas G. Torr 
University of Alabama in Huntsville 
Department of Physics 
Research Institute C-10 
Huntsville, Alabama 35899 


Abstract 

We report the development of a new coating design for applications in the vacuum 
ultraviolet which yields significantly higher reflectivity over selectable bandwidths. We 
demonstrate that the concept can be used to fabricate high performance narrow and 
broadband reflection filters, whose spectral properties can be greatly enhanced by utilizing 
several of these filters in tandem. For example, we have fabricated a narrowband filter 
at the location of the 01 135.6 nm line with 3.2 nm bandwidth, peak transmittance of 
39.3%, and blocking of out-of-band wavelengths is better than 10~ 4 %. The principle of 
our design approach is to use a combination of high (H) - low (L) refractive index dielectric 
pairs such that H + L = A r /2 where H/L < 1. H and L designate the optical thicknesses 
of high and low index film materials. This kind of choice for the high-low ratio reduces 
the effects of absorption for the H films for which the extinction coefficient in the VUV 
is much higher than for the low index film material MgF 2 . The reduced absorption of 
multilayers with H/L < 1, results in a significant increase of reflectivity compared to the 
classical quarterwave stack for which H/L = 1. 
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I. INTRODUCTION 


Narrowband filters that are commercially available in the vacuum ultraviolet wave- 
length region from 120 nm to 160 nm have typical transmittance lower than 15%, and full 
width measured at half of transmittance maximum (FWHM) greater than 25 nm. The 
peak transmittance of the filters centered at the longer wavelengths from 160 nm - 230 
nm lie between 20% and 25% with FWHM > 20 nm. 1 

Malherbe 2 reported the design and the spectral performance of a narrowband filter 
centered at the hydrogen Lyman-a line (121.6 nm) with peak transmittance close to 
15% and FWHM = 9 nm. Blocking of the wavelengths longer than 160 nm is better than 
10 _3 %. The filter has relatively high transmittance for the wavelength region from 126 nm 
- 135 nm; close to 7% at 126 nm and almost 1% at 135 nm. This pass window makes the 
filter not very useful for terestrial imaging applications if the spectral discrimination of the 
neighboring atomic oxygen lines at 130.4 nm and 135.6 nm is desired. A narrowband filter 
centered at 202.5 nm is reported by the same author 3 . The filter has peak transmittance 
greater than 85% and FWHM = 2.5 nm. However, the blocking zone of the filter is very 
short and transmittance for wavelengths longer than 220 nm becomes greater than 85%. 

The calculated and experimental spectral performance of a Fabry-Perot-type nar- 
rowband filter centered at 179 nm was reported by Spiller 4 . The theoretical calculation 
predicted a narrowband filter with resolution A 0 /AA — 60 and peak transmittance of 25%, 
but the measured performance had almost four times smaller resolution and much smaller 
peak transmittance. Discrepancies between the theoretically predicted and experimentally 
obtained spectral curves have been ascribed to the excitation of a surface plasma wave 
traveling along the surface of an aluminum film. 
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A variable bandwidth transmission filter reported by Elias 5 et al had bandwidths 
from 7 nm to 20 nm with peak transmittance from 20% to 40%, respectively. The filter 
was centered at 176 nm, and as in the case of other all-dielectric filters, suffered from pass 
windows in the longer wavelength region. 

Narrowband filters for the vacuum ultraviolet wavelength range (VUV) from 120 nm 
to 230 nm with similar optical properties to those listed above were also reported by 
some other authors 0-9 . The FWHM bandwidths of some filters that we have developed 
recently 10 are smaller than 5 nm with peak transmittance greater than 25% centered at 
wavelengths as low as 135 nm. Since those were all-dielectric filters, the transmittance in 
the longer wavelength region provided insufficient blocking for some applications. 

Currently available VUV broadband filters with bandwidths greater than 10 nm have 
relatively low transmittance and the shape of the transmittance curve is similar to that 
of Fabry-Perot-type filters 1,4-8,11 . For most applications a more rectangular shape for the 
pass band and higher transmittance is required. 

Taking all this into account, it is clear that the VUV spectral range lacks high quality 
narrow and broadband filters such as are available in the visible and infrared parts of the 
spectrum. The reasons for this are certainly a lack of low absorbing film materials in the 
VUV for all-dielectric filters, and the coupling of the incident light into plasma surface 
waves of the metal for metal-dielectric filters. 

Hunter 12 achieved a measure of success in solving this problem by combining two or 
more reflectors in series in order to achieve desired spectral performance for the design 
of vacuum ultraviolet reflection polarizers and analyzers. The idea is that if sufficiently 
high reflectivity can be achieved within the passband, then the in-band exponential loss of 
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reflectivity with additional reflective surfaces becomes insignificant compared to the net 
out-of-band exponential reflectance reduction. However, for the approach to be viable, 
the ratio of in-band to out-of-band reflectivity at each surface should be the order of 10, 
eg. 90% and 10%, respectively. Two reflections, for example, then reduce the in-band 
reflectivity to 81%, whereas the out-of-band reflectivity is reduced to 1%, and so on. 

In this paper we report the successful design and fabrication of narrow and broad- 
band reflection filters with the desired 90% reflectance. These filters are then combined 
into a multiple reflector to provide excellent blocking for out-of-band wavelengths and the 
desired spectral shape for both narrowband and the broadband applications. For exam- 
ple, narrowband filtering with a bandwidth of 3.2 nm and a throughput at the central 
wavelength (135.6 nm) of 39.3% was achieved. Blocking of out-of-band wavelengths is 
better than 10~ 4 %. A multiple reflector was designed to illustrate the performance of a 
broadband filter centered at 175 nm. In this case a bandwidth of 11 nm and throughput 
greater than 54% was realized. Out-of-band blocking is better than 3xl0~ 3 %. 

The next section of this paper reviews the theory of quarterwave (QW), and other 
multilayer options which we have called II filters 13 . The design and spectral performance 
of several narrowband and broadband filters are presented in Section III. Section IV deals 
with experimental techniques including thin film deposition, and substrate handling and 
preparation. A summary is given in Section V, 
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II. ABSORBING MULTILAYERS 


A. Reflectance Transmittance and Absorptance 


The intensity reflection and transmission coefficients for a plane electromagnetic wave 
incident on a multilayer stack are given by 14 


{Mu ± M 13 t7s)t?o - (Mm + MiiVs) 

(M u + MuVs)vo + (M 2 i + M 22 tis)’ 


t = 


2 Vo 


( 2 ) 


(Mu + Mi2rjs)Vo + (Mji + MjjVs)' 
where % and r/s, which are the effective optical functions of the incident medium and the 
substrate, are defined as 


770 = 7lQ COS 0 q, 


(3) 


775 = 715 cos 0S 


(4) 


for s polarization, and as 


COS 00 

Vo = 1 

7lo 


IS 


cos 8 s 

> 

^5 


(5) 

( 6 ) 


for p polarization. Angles, 0 O and 0,, for light propagation through the incident medium 
and the substrate are measured relative to the normal to the film plane, n 0 and n s are the 
refractive indices of the incident medium and the substrate, respectively. It is assumed 
that both the substrate and incident medium have negligible extinction coefficients and 
therefore real optical functions. The terms M tJ ,i, j =1,2 are the elements of the multilayer 
characteristic matrix M which is defined as the product of the matrices of the individual 
layers Mr, where 


( cos Si 
t? 7 r sin Si 


— sin Si \ 
vt 1 ] 

cos Si ) 


(7) 
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The phase thicknesses of the films b\ are given by 


St = ^ Nidi cos ( 8 ) 

*0 

where A 0 is the vacuum wavelength of the incident light. N t is the optical function (usually 
referred to as “optical constant”) of the l-th layer defined as 

Ni = 71/(1 + tKj) = 71| + iniKi = 71/ + ik h (9) 

with K t = k t /n h and where n f is the refractive index, &/ is the extinction coefficient, is 
the physical thickness, and 0/ is the complex angle of the light within the l-th film. 

The reflection and transmission coefficients r and t are complex numbers of the form 

r = |r|e'^ p , (10) 

< = ( 11 ) 

where <f> r and <j> t are the phase changes on reflection and transmission. The intensity 
reflectance i?, transmittance T, and absorptance A of a multilayer are given by 


R = rr*, 

(12) 

T = — «*, 

(13) 

Vo 

A = l-{R + T). 

(14) 


B. Quarterwave Tuned Multilayers 

Multilayer stacks formed by high and low index materials alternating throughout a 
stack are usually referred to as tuned filters or tuned multilayers. The basic design of 
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such a multilayer stack is given in either symmetric [(HL) P H] or asymmetric form [(HL) p j. 
Symmetric [(LH) P L] and asymmetric [(LH) P ] tuned multilayers represent another form 
of the basic design. H and L designate optical thicknesses of high and low refractive 
index film materials respectively, and p is the number of (HL) pairs. The form and the 
symmetry of a stack are determined by the refractive indices of the substrate and incident 
medium 15 . 


If the optical thickness of a film (defined as a product of a physical thickness of a film 
and its refractive index) is equal to one quarter of some reference wavelength A r then this 
film is referred to as a quarterwave (QW) . Applications of QW multilayers are based 
on the fact that for an incident wavelength A 0 = A r the beams reflected from the various 
interfaces will all be in phase, so the reflectance obtained is a maximum. 


We have shown that at zero angle of incidence, the maximum reflectance of a QW 
stack with absorbing film materials is achieved when the number of (HL) pairs p satisfies 
the following condition 10 


P = 




1 + KffKL 



( 15 ) 


where 



nH 


( 16 ) 


and 



( 17 ) 


u fj and tii are refractive indices of high and low index film materials, and k h and k L are 
corresponding extinction coefficients. Values of p > p 0 correspond to numbers of (HL) 
pairs for which both the absorptance and reflectance of the stack are constant and sum 
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to unity i.e. 


R + A->1, (18) 

resulting in a value for the stack transmittance essentially equal to zero. 

An alternative description of the optical properties of a QW stack with absorbing 
films might be given in terms of the standing wave ratio (SWR) 10 . The SWR in general 
is the ratio of the maximum and minimum amplitudes in a standing wave. In multilayer 
theory it is the ratio of the maximum to minimum electric field amplitude in the standing 
wave formed by the interference of the incident and reflected electromagnetic waves. For 
a symmetrical QW stack 17 

SWR = 1 + ^ (19) 

i - Sr 

where R is the intensity reflectance of the stack. The SWR saturated level and therefore 
the maximum reflectance of a QW stack of the form [(HL) P H] is given by 1 * 

Rk = 1 - 2xn 0 — y , (20) 

n H” n L 

where i?#, the ultimate reflectance, is usually referred to as the Koppelmann limit. Thus, 
the reflectance of a QW stack R — ► Rk as p — ► oo, if film materials forming the stack are 
absorbing. For non-absorbing film materials R — ► 1 as p—> oo. 

It should be emphasized that Eq.(20) was derived with some approximations 18 and it 
cannot replace an exact calculation of the maximum reflectance of a QW stack. Equation 
(20) provides an extremely simple relationship between the maximum reflectance obtain- 
able with an absorbing QW stack and the optical properties of film materials forming 
the stack. However, the value obtained for the maximum reflectance must be treated 
as an estimate. In addition, Rk is derived for zero angle of incidence and it does not 
represent the ultimate reflectance for 0 o ^ 0°. Therefore, we shall reference the exactly 
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calculated reflectance [Eq. (12)] as the “Koppelmann limit” for a QW stack at any angle 
of incidence. 

The Koppelmann limit and how to extend it have been the subject of a number of 
theoretical and experimental studies 10-21 . Some approaches involve numerical refinement 
of the film thicknesses in order to reduce the electric field amplitude in the high index 
materials close to the side of incidence. Other authors, assuming values of the extinction 
coefficient lower than 10 -3 , derived approximating formulas for the thicknesses of an 
optimum (HL) pair which would provide the maximum reflectance. 

BaF 2 and LaF 3 are the only available VUV high index film materials identified at 
this time 22 . They both have values of the extinction coefficient greater than 10 -2 for 
wavelengths below 200 nm. This means that most of above mentioned and referenced 
techniques for the improvement of the Koppelmann limit are not applicable to the VUV 
wavelength range. Our approach to this problem is to utilize multilayer stacks with the 
optical thickness ratio H/L < 1 instead of standard QW stacks for which H/L = 1. MgF 3 
film material has values of the extinction coefficient of the order of 10 -< for almost entire 
VUV wavelength region 22 , and it is exclusively used as the low index film material for 
all multilayer designs presented in this paper. Since the extinction coefficient of the high 
index materials (BaF 3 and LaF 3 ) is 100 times greater than that of MgFj, the stacks with 
H/L < 1 have lower absorption and therefore higher reflection than stacks with H/L = 1. 

C. Thirdwave Tuned Multilayers 

The definition of the thirdwave (TW) multilayers is not standardized as in the case 
of the quarterwave stacks. It is defined here as a tuned multilayer which has one material 
with optical thickness equal to one third of some reference wavelength, while the optical 
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thickness of the other material is equal to one sixth of a reference wavelength, i.e. 


n H<^H = 


K 
6 ’ 


riLtii — 


K 
3 ' 


Since k H /k L > 100 for presently available high and low index materials in the VUV (LaF 3 , 
BaF 2 , and MgF 2 ), the choice of which material should have the optical thickness of A r /6 
is obvious. 


The two H and L films in a QW multilayer form an (HL) pair with total optical 
thickness equal to A r /2 which corresponds to a total phase thickness of the pair 6 = w. 
Similarly, an (HL) pair of a TW multilayer has an overall optical thickness equal to A r /2 
which again corresponds to a total phase thickness of the pair equal to w. Thus, the 
principle of a high reflectance at a reference wavelength of a TW tuned multilayer is 
similar to that of a QW stack. The concept of the TW multilayer leads to other kinds of 
multilayers which are discussed in the next section. 

D. n Multilayers 

We define a n multilayer as one whose basic (HL) pair has a total optical thickness of 
A r /2. Optical thicknesses of individual H and L films forming a pair satisfy the following 
condition 

h + l = Y ' (21) 

where A r is the reference wavelength of the multilayer. The total phase thickness of the 
pair is equal to k, i.e. S H + S L = n. Thus, quarterwave and thirdwave stacks are the 
special cases of n multilayers. In a quarterwave stack light reflected from all interfaces is 
in phase, while in a thirdwave and other n multilayers light reflected from each (HL) pair 
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is in phase. Obviously, quarterwave stacks with low-absorbing film materials (available 
in the visible and infrared parts of spectrum) provide higher reflectance with fewer layers 
than other II stacks. However, in the VUV where low-absorbing high index film materials 
do not exist, a II multilayer with smaller physical thickness of H relative to L can provide 
lower absorptance and therefore higher reflectance of the stack. 

Figure 1 illustrates how the Koppelmann limit can be easily extended for zero angle 
of incidence. It shows maximum reflectances calculated for the 99-, 55-, and 35-layer n 
stacks plotted against the ratio of optical thicknesses of high and low index film materials 
(H/L). The reference wavelength is 145 nm, with MgF 2 as the high, LaF 3 as the low index 
materials, and fused silica as the substrate. The Koppelmann limit is R K = 90.8%. The 
99-layer filter has the maximum reflectance of 96.5% when the ratio of optical thicknesses 
H/L = 1/4. This ratio corresponds to the optical thicknesses of H = A r /10 and L = 4A r /10. 
The 55-layer has the maximum reflectance of 95.3% for H/L = 1/3 which corresponds to 
the optical thicknesses of H = A r /8 and L = 3A r /8. The ratio H/L = 1/2 provides the 
maximum reflectance of 93.6% for the 35-layer n stack. From Fig. 1, it follows that at 
zero angle of incidence, stacks with the ratio H/L < 1 provide a significant improvement 
in the maximum reflectance over what can be achieved with QW stacks for which H/L = 
1 . 


Shown in Figures 2 and 3 are maximum reflectances and bandwidths (full width 
measured at half of reflectance maximum) respectively, calculated at 0 O = 45° for the 99-, 
55-, and 35-layer II stacks, plotted against the ratio of optical thicknesses of high and 
low index material (H/L). The reference wavelength is 135.6 nm, with MgF 2 as the low, 
LaF 3 as the high index material, and fused silica as the substrate. The Koppelmann limit 
for the QW stack is R K — 89.5% 0 o = 45° at angle of incidence. The 99-and 55-layer II 
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multilayers (Fig. 2) have a maximum reflectance of 93.0% and 92.7% respectively, when 
the ratio of optical thicknesses is H/L = 1/3. This corresponds to optical thicknesses of 
high and low index film materials 

h- k 

h ~T' 

_ 3A r 
8 

The 35-layer n stack (Fig. 1) has a maximum reflectance of 91.4% for H/L = 1/2. The 
optical thicknesses for this H to L ratio of the n multilayer are given by 



L = 


K 

3 ' 


The lowest decrease of the maximum reflectance occurs for the QW stack when the 
number of films is decreased from 99 to 35. This result is certainly expected since the 
saturation level of the SWR decreases with an increase of the multilayer absorption. Of all 
the multilayers the physical thickness of the high index material (LaF 3 in this particular 
case) within the multilayer is largest for the QW stack. Therefore, the SWR saturation 
level for the QW stacks is achieved with fewer layers compared to other II stacks. As 
should be expected the bandwidth decreases with a decrease of of the ratio H/L. It is 17 
nm for the QW stack and 6.9 nm for the stack with H/L = 1/5 (Fig. 3). 


From Figs. 2 and 3, it follows that for the design of a narrowband reflector centered at 
135.6 nm, the n stack with H/L = 1/4 seems to be the most feasible choice. The 35-layer 
stack with this ratio has optical thicknesses given by 


nHdH = W 


, 4A 

n l d L = ^ 
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Calculated values of the maximum reflectance R max = 88.4% at 135.6 nm and FWHM 
= 7.2 nm. Figure 4 shows the measured and calculated reflectance of this stack with 
fused silica as the substrate material. The measured values of the peak reflectance at 
135.6 nm is 88.3% while the theory predicts the reflectance peak of 88.4%. The measured 
bandwidth of 7.2 nm agrees with the theoretical predictions. 


Figure 5 shows the measured and calculated reflectance at 45° angle of incidence for 
a 29-layer n stack with BaF 2 as the high index material and fused silica as the substrate. 
The ratio H/L = 1/3, and optical thicknesses are given by 


and# = 


A, 

8 


re Z-^L = 


3A r 

8 


The measured reflectance maximum at A r = 135.6 nm is 86.0% with bandwidth of 7.2 
nm. The theory predicts 86.9% maximum reflectance and 6.8 nm bandwidth. 


E. Higher Order n Multilayer Stacks 


The first order QW stacks have wider high reflectance zone than other n multilayers. 
The width of the high reflection zone (AA) WR of a QW multilayer with non-absorbing 
film materials is given by 23 


(AA)tf.K. 


1 

2 (m - 1) + 1 



(24) 


where m is the order of n multilayer. For absorbing film materials the width of the high 
reflection zone is smaller than that calculated using Equation (24). Thus, Equation (24) 
gives a maximum width for a high reflectance zone obtainable for materials with refractive 
indices uh and n^. Thus, an alternative approach for the design of the narrowband 
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reflectors is to utilize the QW or other II stacks with higher order m. However, an increase 
of the order of a n stack from 1 to, say m, changes a total phase thickness of an (HL) pair 
from 7 r — > m7r. The total optical thickness of an (HL) pair is increased by A r /2 when the 
order of a n stack is increased by unity. Increased physical thickness of both high and 
low index film materials results in increased absorptance and therefore lower reflectance 
of the stack. Figure 6 shows the calculated and experimentally obtained reflectance of 
the second order QW stack centered at 135.6 nin at a 45° angle of incidence. BaFj is 
used as the high index material and fused silica as the substrate. The calculated value 
of the maximum reflectance is 83.1% while the measured value is 69.9%. Discrepancies 
between the theory and experiment may be explained by the presence of physical factors 
neglected in the theoretical calculation such as the surface and volume scattering, film and 
substrate contaminations, and film inhomogeneity. The agreement between the theoretical 
prediction and experimentally obtained results is much better for the first order filters 
(Figs. 4 and 5) than for the second order QW stack. This may indicate that thicker films 
in the QW stack are more affected by the physical factors neglected in the multilayer 
theoretical calculation. 
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III. MULTIPLE REFLECTION FILTERS 


A. Narrowband 

Reflection filters shown in Figures 4 and 5 have measured reflectances at the central 
wavelength of 86.0% and 88.3%, respectively. They both have relatively high reflectance in 
the shorter wavelength region and an average reflectance of the order of 10% for the longer 
wavelengths. The reflectance outside the filter’s “pass” zone and the bandwidth of the pass 
zone can be reduced by means of multiple reflections from two or more reflection filters 
upon which light is incident at an angle of 45°. The overall transmissions of combinations 
of four and six 29-layer filters are shown in Figures 7 and 8, respectively. 

The maximum throughput of the 4-filter combination at the central wavelength A = 
135.6 nm is 53.7% and the bandwidth of 4.3 nm. The blocking for shorter wavelengths is 
better than 0.7% and for longer wavelengths is better than 10~ 2 %. If further improvement 
of the blocking for wavelengths outside the pass zone is required, then 6 or more filters 
can be combined. However, adding more filters will reduce the overall transmittance of 
the combination. The maximum throughput of the 6-filter combination at the central 
wavelength is 39.3% with the bandwidth of 3.2 nm. The blocking outside the pass zone is 
better than 10 -4 % for longer and better than 0.07% for shorter wavelengths. If a 4 mm 
thick MgF 2 window is placed at the entrance of the combination than the transmittance 
of shorter wavelengths may be less than 10" 4 %. 

B. Broadband 

The pass zone of a broadband filter is bounded by a lower and upper wavelengths. 
Ideally, the spectral components of the incident light, with wavelengths shorter than the 
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lower and longer than the upper wavelength of the filter, together referred to as the out- 
of-band spectrum, are rejected. In the design examples that follow, wavelengths of the 
out-of-band spectrum are rejected by means of multiple reflections from QW stacks. The 
rejection at shorter wavelengths might be improved by the suitable choice of the window 
material, placed at the entrance of a multireflector combination. Windows made of BaF 3 
and CaF 2 absorb wavelengths below 135 nm, and 125 nm respectively, while fused silica 
and AI2O3 may be used for broadband filters with a lower pass limit above 145 nm 13,34 . 

Because of the narrower high reflection zone the n multilayers with H/L < 1 are not 
suitable for the design of broadband reflectors. QW stacks for which H/L = 1 are a 
better choice for the design and fabrication of these filters. Figure 9 shows the measured 
and calculated reflectance at 45° incident angle for the 25-layer QW stack. The reference 
wavelength is 175 nm, LaF 3 is used as the high index material, and fused silica as the sub- 
strate. The measured peak value of the reflectance is 90.3% and the bandwidth measured 
at half of the reflectance maximum is 19 nm. The average reflectance for the out-of-band 
wavelengths is of the order of 10%. The Koppelmann limit for LaF 3 and MgF 2 at this 
wavelength is R K = 94.4% , and the theoretically predicted value of the peak reflectance 
is 91.6%. 

The agreement between the theory and the experiment is much better than for the 
second order QW stack (Fig. 6). Again this may indicate that thicker films in the second 
order QW stack are more affected by the physical factors neglected in the multilayer 
theoretical calculation. 

Figure 10 shows the transmittance of combinations of four and six of these reflectors. 
The 4-reflector combination has a peak transmittance of 66.3% and the bandwidth of 
11.5 nm. The rejection of the out-of-band spectrum up to 300 nm is better than 0.1%. 
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If a 4 mm thick fused silica parallel plate is used as the window then the rejection of the 
shorter wavelengths of the out-of-band spectrum is better than 10~ 3 %. The 6-reflector 
combination has a peak transmittance of 54.04% with bandwidth of 11 nm. This com- 
bination provides much better rejection for out-of-band wavelengths. The transmittance 
of the combination for out-of-band wavelengths is less than 3 x 10~ 3 %. If a 4 mm thick 
fused silica is used as the window than the transmittance of the combination for shorter 
wavelengths of the out-of-band spectrum is less than 10 -5 %. 
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IV. THE EXPERIMENTAL ARRANGEMENT 


All depositions are made on 12.7 mm diameter by 2 mm thick MgF 2 and fused silica 
substrates with root mean square roughness less than 1 nm. The substrates are cleaned by 
the supplier (Acton Research Corp., Acton, Massachusetts) using the following procedure: 
optical soap wash, water rinse, ethanol soak then ultrasonic bath, fresh ethanol rinse, and 
finally a Freon rinse. The substrates are packed in a clean room environment in Delrin 
holders. They are removed from Delrin holders immediately prior to mounting in the 
substrate holder for deposition. Mounting in the substrate holder is done in a laminar 
flow bench. During the transport from the bench to the vacuum chamber the substrate 
is kept in a flow of dry nitrogen. All depositions are made at the University of Alabama 
in Huntsville, Optical Aeronomy Laboratory and spectrophotometric measurements are 
made at the Atomic Physics Branch of the NASA/Marshall Space Flight Center. 

The vacuum system consists of a cryo-pump and a sorption pump giving an oil- 
free environment for all depositions and therefore providing a very low probability for 
hydrocarbon contamination of the films. The film materials BaF 2 and LaF 3 are prepared 
for vacuum deposition by CERAC with a typical purity of 99.9% while MgF 2 (99.95%) is 
standard BALZERS coating material. The films are deposited with low deposition rates 
of 0.18 - 0.20 nm/sec on the substrates heated to 200°C. The pre-deposition pressure is 
7.5xl0~ 7 while the deposition pressure varies between 1.7 and 2.3xl0‘ 6 for all coating 
materials. The temperature is monitored with a Chromel-Alumel thermocouple attached 
to an aluminum ring substrate holder. The substrate and the ring holder are placed in 
the center of a 40 mm diameter by 6 mm thick stainless steel plate. A quartz crystal 
is used for the film thickness and deposition rate monitoring. The depositions are made 
with an electron gun. The gun has fixed voltage of 10 KV and low power depositions are 
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maintained by supplying low current to the filament. The source to the substrate distance 
is 500 mm and the source to the monitor distance is 350 mm. 
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V. SUMMARY 


The idea of utilizing multiple reflections from n multilayer reflection filters upon which 
the light is incident at 45° angle constitutes the basis of our narrowband and broadband 
filters. The multiple reflector combinations provide superior spectral performance for both 
the narrowband and the broadband filters over what is currently available. 

The high reflective multilayer coatings for the narrowband reflection filters are de- 
signed using II multilayer stacks with H/L < 1 instead of the classical QW stacks for 
which H/L = 1. It is shown that the II stacks with H/L < 1 have higher reflectance and 
smaller bandwidth than the QW stacks, and are therefore more suitable for the design of 
narrowband reflection filters. It is demonstrated that the QW stack maximum reflection 
limit at both 0 O = 0° (Koppelmann limit ) and 0 o = 45° can be easily extended for coating 
materials available in the VUV. A peak reflectance of more than 96% for zero angle of 
incidence can be achieved for wavelengths as low as 135.6 nm when II stacks with the 
ratio H/L < 1 are utilized. Stacks with more than 99 layers can provide further improve- 
ments of the maximum reflectance, but the lowest value of the H/L ratio is limited by the 
smallest feasible thickness of the high index films. Furthermore, II multilayers with H/L 
< 1 are not suitable for broadband reflectors since they generally do not have a width of 
the high reflectance zone as wide as the QW stacks. 

The agreement between the theoretical and experimental results is much better for 
first order filters (Figs. 4, and 5) than for the second order QW stack (Fig. 6) and for 
QW stacks centered at longer wavelengths (Fig. 9). This may indicate that thicker films 
in the QW stack are more affected by the physical factors neglected in the multilayer 
theoretical calculation such as the surface and volume scattering, film inhomogeneity, and 
film contamination. 
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FIGURE CAPTIONS 


Figure 1: The maximum reflectance of the n stack calculated for zero angle of incidence 
at 145 nm. Diamonds represent 99-layer, triangles 55-layer, and squares 35-layer stacks. 
H = Lanthanum Fluoride, and L = Magnesium Fluoride. The Koppelmann limit is 90.8%. 

Figure 2: The maximum reflectance of the n stack calculated for 45 degrees angle of 

incidence at 135.6 nm. Diamonds represent 99-layer, triangles 55-layer, and squares 35- 
layer stacks. H = Lanthanum Fluoride, and L = Magnesium Fluoride. The Koppelmann 
limit is 89.5%. 

Figure 3: The full width at half of the reflectance maximum of the II stacks calculated 

for 45 degrees angle of incidence. The maximum reflectance of the stacks as a function of 
the H/L ratio is shown in Fig. 2. 

Figure 4: The measured (squares) and calculated (diamonds) reflectance of the 35-layer 

n stack for 45 degrees angle of incidence centered at 135.6 nm. The optical thickness ratio 
H/L = 1/4 with H = Lanthanum Fluoride, and L = Magnesium Fluoride. 

Figure 5: The measured (squares) and calculated (diamonds) reflectance of the 29-layer 

n stack for 45 degrees angle of incidence centered at 135.6 nm. The optical thickness ratio 
H/L = 1/3 with H = Barium Fluoride, and L = Magnesium Fluoride. 

Figure 6: The measured (squares) and calculated (diamonds) reflectance of the 35- 

layer second order QW stack for 45 degrees angle of incidence centered at 135.6 nm. H 
= Barium Fluoride, and L = Magnesium Fluoride. 
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FIGURE CAPTIONS 

(Continued) 


Figure 7: The transmittance of the combination of four 29-layer filters shown in Figure 

5. The bandwidth is 4.3 nm, and a peak transmittance at 135.6 nm is 53.7%. 

Figure 8: The transmittance of the combination of six 29-layer filters shown in Figure 

5. The bandwidth is 3.2 nm, and a peak transmittance at 135.6 nm is 39.3%. 

Figure 9: The measured (squares) and calculated (diamonds) reflectance of the 25-layer 
QW stack for 45 degrees angle of incidence centered at 175 nm. H = Lanthanum Fluoride, 
and L = Magnesium Fluoride. 

Figure 10: The transmittance of combinations of four (solid line) and six (dashed line) 

25-layer QW stacks shown in Fig. 9. 
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degrees angle of incidence at 135.6 nm. Diamonds represent 99-layer, 
triangles 55-layer, and squares 35-layer stacks. H = Lanthanum Fluoride, 
and L = Magnesium Fluoride. The Koppelmann limit is 89 5%. 
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Figure 3: The full width at half of the reflectance maximum of the n 

stacks calculated for 45 degrees angle of incidence. The maximum reflectance 
of the stacks as a function of the H/L ratio is shown in Fig. 2. 
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Figure 4: The measured (squares) and calculated (diamonds) reflectance 
of the 35-layer 11 stack for 45 degrees angle of incidence centered at 
135.6 nm. The optical thickness ratio H/L = 0 25 with H = Lanthanum 
Fluoride, and L = Magnesium Fluoride. 
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135.6 nm. The optical thickness ratio H/L = 1/3 with H = Barium Fluoride, 
and L = Magnesium Fluoride. 
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Figure 6: The maesured (squares) and calculated (diamonds) reflectance 
the 35— layer second order QW stack for 45 degrees angle of incidence 
centered at 135.6 nm. H = Barium Fluoride and L = Magnesium Fluoride 
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Figure 8: The transmittance of the combination of six 29-layer filters 
shown in Fig. 5. The bandwidth is 3.2 nm, and a peak transmittance at 
135.6 nm is 39.3%. 
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Figure 9: The measured (squares) and calculated (diamonds) reflectance 
of the 25-layer QW stack for 45 degrees angle of incidence centered at 
175 nm. H = Lanthanum Fluoride, and L = Magnesium Fluoride. 
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Figure 10: The transmittance of combinations of four (solid line) 
and six (dashed line) 25-layer QW stacks shown in Fig. 9. 



APPENDIX C 



High Throughput Narrowband 83.4 nm Self-Filtering Camera 


Muamer Zukic and Douglas G. Ton- 
University of Alabama in Huntsville 
Research Institute C-10 
Huntsville, AL 35899 

Phone: (205) 895-6238 FAX: (205) 895-6717 
and 

Marsha R. Torr and James F. Spann 
NASA/Marshall Space Flight Center 
Huntsville, AL 35812 

Phone: (205) 544-7591 FAX: (205) 544-5862 
Abstract 

Photometric imaging of ionospheric/magnetospheric 0 II emission at 83.4 nm is a 
primary objective for imaging and tomographic mapping of magnetospheric 0 + ions. 
However, instrumental sensitivity has been a major barrier to realizing this goal. We 
report an instrumental design employing a low focal ratio (f/# less than 2) three- 
mirror camera where the reflecting surfaces act as both narrow-band reflection filters 
at 83.4 nm and as a high quality imaging system. The design includes multilayer 
coatings with reflectances that are relatively insensitive to the angle of incident light (0 
50°). The peak reflectance per mirror is more than 60% at 83.4 nm with the average 
reflectance for out-of-band wavelengths of less than 5%. The net reflective transmission 
for the three mirrors is greater than 20% with 5nm bandwidth and 10“ 2 % rejection for 
out-of-band wavelengths increasing to 5x10-“% at H Ly a (121.6 nm). The low focal 
ratio combined with high system transmission results in about 2 orders of magnitude 
improvement in sensitivity over previously reported designs. The camera spot size 
throughout the entire image plane is less than 50 /mi for an 8° field of view, improving 
with smaller fields of view. 
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